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ABSTRACT 


Filmwise  condensation  heat-transfer  measurements  of 
steam  were  made  on  horizontal  tubes  under  vacuum  and  near- 
atmospheric  pressures.   Data  were  taken  for  a  smooth  tube 
and  for  21  tubes  which  contained  rectangularly-shaped,  low 
integral  fins.   The  fin  geometry  was  systematically  varied 
in  order  to  investigate  the  dependence  of  the  steam-side 
heat-transfer  coefficient  on  fin  spacing,  thickness  and 
height. 

The  condensation  process  was  found  to  be  most  sensitive 
to  fin  spacing  and  to  be  weakly  dependent  on  fin  thickness. 
When  the  fin  height  was  increased  from  1.0  to  2.0  mm,  the 
increase  in  the  heat-transfer  coefficient  was  about  20 
percent,  while  the  increase  in  condensing  area  was  51 
percent. 

The  best  performing  finned  tube  had  a  fin  spacing  of  1.5 
mm,  a  fin  thickness  of  1.0  mm  and  a  fin  height  of  2.0  mm. 
It  showed  an  enhancement  over  the  smooth  tube  of  between  4 
to  5  under  vacuum  conditions  and  around  6  at  atmospheric 
conditions.   Thus,  the  use  of  finned  tubes  may  result  in  a 
significant  reduction  in  the  size  and  weight  of  present-day 
steam  condensers. 
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Bottom  of  the  Tube  to  the  Position  at  Which 
the  condensate  Just  Fills  the  Interfin  Space 
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I.   INTRODUCTION 

Modern  marine  propulsion  plants  must  be  designed  to 
provide  maximum  power  with  small  size,  low  weight  and 
reasonable  cost.   A  heavy  power  plant  may  have  an  impact  on 
ship  dimensions  and  displacement  since  the  power  required  to 
achieve  a  given  speed  increases  with  the  displacement.   A 
sizeable  propulsion  plant  will  occupy  space  that  is  valua- 
ble, while  the  machinery  room  will  have  poor  accessibility. 

A  significant  decrease  in  size  and  weight  of  a  steam 
power  plant  can  be  achieved  by  reducing  the  size  of  the 
condenser.   Present  condensers  are  large  because  of  the  use 
of  smooth  tubes  and  generous  design  margins.   The  heat- 
transfer  performance  of  any  condenser  tube  is  limited  by  the 
thermal  resistances  on  the  water-side,  in  the  wall  and  on 
the  steam-side  of  the  tube.   Reduction  of  any  of  these 
resistances  will  contribute  to  a  lighter  and  smaller 
condenser . 

In  recent  years,  various  efforts  have  been  made  in  order 
to  increase  the  outside  (i.e.,  steam-side)  heat-transfer 
coefficient,  and  these  methods  have  been  thoroughly  reviewed 
by  Marto  [Ref.  1] .   Outside  enhancement  techniques  include 
the  use  of  low  integral  fins,  roped  tubes,  fluted  tubes, 
drainage  strips  attached  on  the  tubes  and  applied  coatings 
to  promote  dropwise  condensation. 
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A.   FILMWISE  CONDENSATION  ON  FINNED  TUBES 

Filmwise  condensation  heat  transfer  on  a  smooth  tube  is 
limited  by  the  thickness  of  the  condensate  layer.   This 
layer  forms  a  relatively  high  conduction  resistance  between 
the  vapor  and  the  tube  surface.   To  enhance  filmwise  conden- 
sation heat  transfer,  it  is  then  important  to  thin  this 
condensate  film.   One  common  way  to  decrease  this  condensate 
film  is  to  use  radial  fins  attached  to  the  external  surface 
of  the  tube  in  order  to  promote  surface-tension  effects. 

The  theoretical  treatment  of  this  problem  is  very  diffi- 
cult owing  to  the  large  number  of  controlling  parameters 
such  as  the  three-dimensional  condensate  flow  pattern,  vapor 
shear,  surface-tension  forces,  wall  conduction  effects  and 
condensate  retention.   Studies  following  Gregorig  [Ref.  2] 
have  demonstrated  the  important  role  of  surface  tension  in 
modifying  the  condensate  film  over  the  fin  surface.   The 
presence  of  fins  causes  three  main  effects: 

1.  the  condensing  area  increases, 

2.  the  surface-tension  forces  thin  the  condensate  film 
over  the  tip  of  the  fins  causing  improved  heat 
transfer,  and 

3.  the  fins  may  cause  condensate  retention  especially 
on  the  lower  part  of  the  tube. 

The  first  two  effects  are  beneficial,  while  the  third 
undermines  the  heat-transfer  performance.   The  second  effect 
is  most  important  in  the  upper  portion  of  the  tube,  where 
there  is  little  condensate  retention.   For  steam  condensa- 
tion data  on  finned  tubes,  Yau  et  al.  [Ref.  3]  measured 
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enhancement  ratios  (i.e.,  the  ratio  of  vapor-side  coeffi- 
cient of  a  finned  tube  to  the  value  of  smooth  tube)  in 
excess  of  the  increase  in  condensing  area.   This  additional 
enhancement  is  most  probably  due  to  the  thinning  effect  of 
the  surface-tension  forces  acting  on  the  condensate  film. 
This  phenomenon  can  be  explained  by  referring  to  Figure  1.1 


Fin 


Condensate 


Figure  1.1.   Condensate  Profile  at  Top  Part  of  Tube 

Since  the  condensate  film  has  a  concave  shape  at  point 
A,  the  pressure  within  the  condensate  at  this  point  is 
smaller  than  the  local  vapor  pressure,  and  is  given  by: 


Pv  "  T, 


(1.1) 


A  similar  expression  may  be  written  at  point  C 


=   P 


a 


v 


(1.2) 
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Note  that  the  local  radii  of  curvature  at  points  A  and  C  are 
very  large;  therefore,  the  liquid  pressure  at  these  points 
nearly  equals  the  vapor  pressure.   On  the  other  hand,  at 
point  B  the  condensate  profile  has  a  more  concave  shape, 
thus  resulting  in  a  lower  pressure  than  at  either  points  A 
or  C: 


PB   =   Pv  -■£-  (1.3) 

B 


As  a  result,  there  exists  a  pressure  difference  from  point  A 
to  point  B  as  given  by: 


AP_n   =   a(— -)  (1.4) 

**        rB    rA 


and  in  a  similar  way  from  point  C  to  point  B  as  given  by: 


APCB   "   °<^-  ^>  (1-5» 


In  general,  these  pressure  differences  are  positive,  causing 
the  condensate  to  be  pulled  toward  point  B  (which  acts  as  a 
condensate  run-off  channel) ,  thus  thinning  the  condensate 
film  on  the  fins  as  well  as  in  the  interfin  area.   This 
thinning  effect,  therefore,  increases  the  heat-transfer 
performance  in  excess  of  the  surface  area  gain. 
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On  the  other  hand,  as  mentioned  above,  the  third  effect 
of  the  fins  which  is  to  retain  the  condensate,  especially  on 
the  lower  part  of  the  tube,  creates  a  relatively  thick 
condensate  film,  causing  lower  heat  transfer.   When  the  fin 
spacing  is  reduced,  a  high  surface-tension  condensate  will 
bridge  the  fins;  therefore,  the  heat-transfer  performance 
may  decrease  despite  that  the  condensing  area  has  increased. 
Thus,  owing  to  the  competing  nature  of  these  effects,  the 
selection  of  a  best  fin  geometry  needs  considerable  attention. 

The  effect  of  fin  spacing  and  fluid  properties  (such  as 
surface  tension)  on  condensate  retention  has  been  well 
studied  by  various  investigators  [Refs.  4,5,6].   A  detailed 
discussion  of  these  studies  is  presented  in  the  next  section. 

One  more  important  effect  that  is  common  for  both  smooth 
and  finned  tubes  is  the  oscillation  of  the  condensate  layer 
resulting  from  the  intermittent  departure  of  condensate 
droplets  from  the  lower  part  of  the  tube  surface.   When  a 
droplet  departs,  the  nearby  film  stretches  (becomes  thinner) 
and  higher  heat  transfer  takes  place  locally.   After  the 
droplet  has  departed,  the  film  grows  again  and  becomes 
thicker.   The  frequency  and  amplitude  of  this  oscillation 
are  functions  of  the  heat  flux,  fin  geometry  and  the  proper- 
ties of  the  fluid.   To  the  best  of  this  author's  knowledge, 
this  transient  phenomenon  has  not  yet  been  investigated. 
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B.   SURVEY  OF  ANALYTICAL  MODELS  AND  EXPERIMENTAL  WORK 

In  1948,  Beatty  and  Katz  [Ref.  7]  developed  a  simple 
analytical  model  to  predict  the  condensing  heat-transfer 
coefficient  on  a  horizontal  finned  tube.   Their  model  was 
supported  by  various  other  researchers,  as  stated  in 
[Ref.  8].   The  final  form  of  the  Beatty  and  Katz  prediction 
has  the  following  equations: 

(kfP  (p  -p  )gh   )  1/4     1/4 
hBk  =   0.689  J_       fg      [A]  (1.6) 

(        yfAT      )      e 


where 


1/4       A  4>A 

[tf]  =    a^~  ~rn  +  1-3  a— L  -T74  '  (1'7) 

e  eff  D1/q  Aeff  X1/4 

o 


A    =   7T  D  L  (1.8) 

r        o 


Af   =   |(D^  -  D^)NL  ,  (1.9) 


A  «   =   A   +  d>A._#  and  (1.10) 

eff      r     f 

2     2 

tt(D;  -  D^) 

x  =  f4D     •  (LID 


The  major  fault  with  their  model  was  the  neglect  of  any 
surface-tension  effects.   They  therefore  assumed  that  the 


22 


Nusselt  equation  can  be  applied  to  all  the  external  surfaces 
of  the  tube  (i.e.,  tip,  sides  of  fins  and  interfin  space). 
According  to  [Ref.  8  and  9],  this  model  overpredicts  the 
condensing  coefficient  especially  for  those  fluids  that  have 
high  surface  tensions  (for  example,  water). 

Rudy  et  al.  [Ref.  5]  attempted  to  include  the  surface- 
tension  effect  by  modifying  the  original  Beatty  and  Katz 
equation.   To  obtain  this  modification,  they  first  investi- 
gated the  retention  angle  problem.   They  performed  static 
(i.e.,  with  no  condensation)  and  dynamic  (i.e.,  with  conden- 
sation) measurements  of  the  retention  angle  on  tubes  with 
fin  densities  from  0.75  to  1.38  fins/mm.   They  found  that 
the  retention  angle  (see  Figure  1.2)  can  be  very  well  pre- 
dicted (within  ±10  percent)  by  the  following  equation: 


*   =   cos"1[1  "  j£kq]  (1-12» 


The  above  equation  was  tested  against  many  fluids  such  as 
ref rigerant-11,  pentane  and  water.   It  is  worth  mentioning 
that  Honda  et  al.  [Ref.  6]  and  Owen  et  al.  [Ref.  8], 
obtained  essentially  the  same  equation  by  using  different 
models.   Next,  Rudy  et  al.  modelled  the  condensation  process 
assuming  no  heat  transfer  in  the  flooded  portion  of  the 
tube.   They  developed  the  following  prediction: 


hRW   "   (1  "  ?>hBK  (1-13) 
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External  Diameter  of  fins 


Root  Diameter  of  fins 


Tube  Wall 


Retained  Condensate 


Rentention  Angle  (\jj) 


Figure  1.2.   Typical  Schematic  of  Retention  Angle 
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where  h    is  given  by  equation  (1.6).   Unfortunately,  this 
model  tends  to  underpredict  most  of  the  available  experimental 
data  by  as  much  as  about  60  percent  [Ref.  9J. 

Owen  et  al.  [Ref.  8]  modified  the  Rudy  et  al.  model  by 
including  heat  transfer  from  the  flooded  part  of  the  tube. 
For  the  fluids  that  they  tested,  their  model  predicted  the 
condensing  coefficient  within  ±30  percent.  The  prediction 
that  they  developed  is  based  on  the  following  equations : 


h    =   (1  -  ^)h   +  ^hT  (1.14) 

OW  IT    U     TT  L 


725  ( ^f^tl^v^fa ) 1/4  ^^  . 

{     yfAf      ) 


h   =   0.725  !-±-i i v- ^      l^)  ;  (1.15) 

U  I  "T^T  (  D 


hL   ■    (irr-  it'"1  (1-16> 

ef  f     c 


where 


h    -   0.725ik^(^-Pv)gh^j1/4[A11/4, 

j        UfAT        )       Df 


and    (1.17) 


heff   =  ^r-  ■  (1-18) 


where 


k  c-     =   (l-sN)k-.   +  sNk^  (1.19) 

eff  fin      f 
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The  above  analysis  uses  Beatty  and  Katz  prediction  (equation 
(1.6))  for  the  unflooded  part  of  the  tube.   As  mentioned 
earlier,  the  surface-tension  effect  is  not  included  in  equa- 
tion (1.6).   However,  as  discussed  in  the  previous  section 
(Section  A) ,  the  surface-tension  forces  play  an  important 
role  in  filmwise  condensation  by  thinning  the  condensate 
film.   Therefore,  the  Owen  model  needs  further  improvements 
to  account  for  the  surface-tension  effect  at  the  unflooded 
part  of  the  tube. 

Recently,  Honda  et  al.  [Ref.  9],  developed  a  model  to 
predict  the  condensing  coefficient  by  accounting  for  surface- 
tension  effects  and  for  the  non-isothermal  behavior  of  the 
fins.   Their  model  can  predict  most  of  the  available  experi- 
mental data  including  11  fluids  and  22  tubes  within  ±20 
percent.   Honda  et  al.  provide  a  comparison  of  these  four 
models  with  experimental  data.   It  is  worth  mentioning  that 
the  worst  disagreement  (within  ±30  percent)  between  the  data 
and  their  model  occurs  for  the  steam  data  of  Rose  [Ref.  3]. 
This  is  probably  owing  to  the  high  surface  tension  of  water, 
as  the  surface-tension  effects  are  still  inadequately 
accounted  for  even  in  this  model.   Further,  it  is  shown  that 
all  the  models  except  that  of  Honda  overpredict  these  steam 
data.   In  comparison  with  the  other  models,  the  Honda  model 
shows  considerable  promise.   However,  it  must  be  noted  that 
this  model  is  rather  laborious  as  it  requires  considerable 
computer  time  for  iterative  finite  difference  techniques. 
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The  Honda  model  may  require  further  improvements  to  predict 
steam  data  more  accurately.   Also,  most  importantly,  a  more 
practical  model  must  be  developed  so  that  the  designers  may 
use  it  with  confidence.   Such  a  model  may  consist  of  a 
combination  of  a  few  correlations.   Thus,  it  is  extremely 
important  to  generate  a  sufficient  amount  of  experimental 
data  covering  all  important  geometric  variables  (such  as  fin 
spacing,  fin  thickness,  fin  height,  fin  shape,  tube  diameter, 
etc.)  for  condensation  of  steam  on  finned  tubes. 

In  1961,  Staub  [Ref.  10]  conducted  steam  condensation 
experiments  with  horizontal  fluted  tubes  under  vacuum  and 
atmospheric  conditions.   He  used  a  bell-shaped  (rather  than 
rectangular)  fin  profile  on  a  root  diameter  of  14.9  mm.   These 
fins  had  a  pitch  of  1.0  mm,  a  thickness  of  0.4  mm  at  the 
root  of  the  fin  and  a  height  of  0.5  mm.   He  reported  a  vapor- 
side  heat-transfer-coefficient  enhancement  (i.e.,  the  ratio  of 
vapor-side  coefficient  of  finned  tube  to  the  value  of  smooth 
tube)  of  around  3  for  both  vacuum  and  atmospheric  pressure 
conditions . 

At  Queen  Mary  College  of  London,  Yau  et  al.  [Ref.  3]  obtained 
steam  condensation  data  at  atmospheric  pressure  on  various 
rectangularly-shaped  finned  tubes  having  a  root  diameter  of 
12.7  mm.   They  found  an  optimum  spacing  of  1.5  mm  for  a  fin 
thickness  of  0.5  mm  and  a  fin  height  of  1  mm.   Their  optimum 
tube  resulted  in  a  vapor-side  heat- transfer  coefficient 
enhancement  of  3  to  4 . 

On-going  research  is  underway  at  the  Naval  Postgraduate 
School  to  study  the  feasibility  of  using  low  integral-finned 
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tubes  with  steam.   A  test  apparatus  has  been  constructed  by 
Krohn  [Ref.  11]  to  systematically  study  steam  condensation 
on  a  single  horizontal  tube.   The  test  apparatus  was  instru- 
mented and  tested  by  Graber  [Ref.  12],  and  further  system 
improvements  for  leak  tightness  were  made  by  Poole  [Ref.  13]. 
Poole  conducted  experiments  both  under  vacuum  and  at  pressures 
near-atmospheric.   He  used  machined,  rectangular  fins  with 
a  fin  thickness  and  a  fin  height  of  1  mm  on  a  root  diameter 
of  19  mm.   He  tested  a  total  of  six  finned  tubes  with  fin 
spacings  of  0.5,  1.0,  1.5,  2.0  ,  4.0  and  9.0  mm  as  well  as  a 
smooth  tube.   His  results  showed  that  the  best  enhancement 
of  3  to  4  occurred  for  the  tube  with  a  fin  spacing  of  1.5  mm. 
However,  Poole  had  problems  with  partial  dropwise  condensa- 
tion since  the  total  filmwise  mode  lasted  only  for  about  20 
minutes  during  his  experimental  runs.   Thus,  it  became 
necessary  to  repeat  his  test  runs  in  order  to  verify  his 
preliminary  results. 

C.   OBJECTIVES 

The  main  objectives  of  this  research  effort  were,  there- 
fore, to: 

1)  establish  a  preparation  procedure  that  will  ensure 
filmwise  condensation  on  copper  finned  tubes. 

2)  take  data  on  an  instrumented  smooth  tube  to  determine 
the  water-side  heat-transfer  coefficient  with  internal 
enhancement. 
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3)  take  data  on  a  smooth  tube  to  check  the  data-reduction 
program  and  to  establish  a  basis  for  comparison  with  the 
finned  tubes. 

4)  take  data  on  tubes  with  rectangularly  shaped  fins  of 
various  fin  spacings,  fin  thicknesses  and  fin  heights  to 
obtain  a  possible  optimum  fin  geometry. 
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II.   DESCRIPTION  OF  APPARATUS 

A.   SYSTEM  OVERVIEW 

The  apparatus  used  for  this  research  was  essentially  the 
same  used  in  [Ref.  13],  with  several  small  modifications.   A 
schematic  of  the  system  is  shown  in  Figure  2.1.   Steam  was 
generated  in  a  304.8  mm  (12  in)  diameter  Pyrex  glass  boiler 
using  ten  4000  Watt,  480  Volt  Watlow  immersion  heaters. 
Passing  through  a  304.8  mm  (12  in)  to  152.4  mm  (6  in) 
reducing  section,  the  steam  travelled  upward  through  a  Pyrex 
glass  section  2.44  m  (8  ft)  in  length,  around  a  180-degree 
bend  and  back  down  a  straightening  section  1.52  m  (5  ft)  in 
length  before  entering  the  stainless-steel  test  section. 
The  condenser  tube  to  be  tested  was  mounted  horizontally  in 
the  test  section  behind  a  viewport  to  permit  visual  observa- 
tion of  the  condensing  process.   Figure  2.2  shows  the 
details  of  the  test  tube  mounted  in  the  test  section. 

Steam  that  did  not  condense  on  the  test  tube  passed  into 
a  stainless  steel  auxiliary  condenser,  and  all  condensate 
was  returned  via  gravity  to  the  boiler.   The  auxiliary 
condenser  was  constructed  of  two  9.5  mm  (3/8  in)  diameter, 
water-cooled  copper  lines  helically  coiled  to  a  height  of 
457  mm  (18  in) . 

Cooling  water  for  the  test  tube  was  provided,  on  a  once- 
through  basis,  by  two  centrifugal  pumps  connected  in  series. 
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Water  to  these  pumps  was  provided  through  a  large  sump 

3 
(capacity  about  0.4  m  ),  which  was  filled  with  a  continuous 

supply  of  filtered  tap  water.   The  water  could  be  throttled 

from  zero  to  0.55  1/s  (8.8  gpm) ,  giving  a  maximum  water 

velocity  of  4.4  m/s  (14.4  ft/sec)  through  the  tube. 

A  continuous  supply  of  tap  water  was  used  for  cooling 
the  auxiliary  condenser.   Throttling  the  flow  of  tap  water 
through  the  auxiliary  condenser  was  the  means  used  to  vary 
the  internal  pressure  of  the  test  apparatus. 

A  vacuum  pump  was  used  to  create  vacuum  and  for  removal 
of  noncondensing  gases  from  the  test  section.   The  vacuum 
pump  was  in  operation  throughout  the  data  collection  time 
interval. 

B.   SYSTEM  INSTRUMENTATION 

The  input  power  to  the  heaters  was  changed  through  a 
panel-mounted  potentiometer.   A  converter  with  an  input 
voltage  of  440  VAC  generated  a  signal,  which  was  fed  to  the 
data  acquisition  system  in  order  to  calculate  the  power 
input  to  the  heaters.   A  complete  explanation  of  the  con- 
verter function  is  given  by  Poole  [Ref.  13]. 

The  internal  pressure  of  the  system  was  measured  manually 
by  a  U-tube,  mercury-in-glass  manometer  graduated  in  milli- 
meters.  Unavoidably,  steam  would  condense  in  the  manometer. 
Therefore,  the  varying  height  of  the  water  column  in  the 
manometer  needed  to  be  accounted  for  when  measuring  the 
system  pressure. 
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All  temperatures  were  measured  using  0.25  mm  diameter, 
copper-constantan  thermocouples.   Two  thermocouples  were  used 
for  the  steam,  one  for  the  condensate  return  and  one  for  the 
ambient.   The  calibration  procedure  for  these  thermocouples 
is  described  in  [Ref .  13] .   The  temperature  rise  of  the  cooling 
water  through  the  test  section  was  measured  by  two  Hewlett- 
Packard  (HP)  2804A  quartz  thermometers  as  well  as  by  a  10- 
junction,  series-connected,  copper-constantan  thermopile. 
The  temperatures  were  fed  directly  into  the  data-acquisition 
system  as  described  below.   During  this  entire  investigation, 
these  two  independent  means  of  cooling-water-temperature  rise 
measurements  agreed  to  within  ±0.02  K.   Note  that  the  coolant 
temperature  rise  was  as  low  as  0.7  K  for  the  smooth  tube 
(under  vacuum  conditions),  while  it  was  in  excess  of  1.2  K 
for  the  finned  tubes.   Thus,  the  coolant  temperature  rise  was 
measured  with  an  uncertainty  better  than  ±3  percent  for  finned 
tubes . 

C.   DATA  ACQUISITION  SYSTEM 

An  HP  3497A  Data  Acquisition  System,  as  controlled  by 
an  HP  9826A  computer,  was  used  to  monitor  the  system  tempera- 
tures and  the  converter's  signal  representing  the  power  input 
to  the  heaters .   The  data  were  immediately  processed  and  also 
were  stored  on  computer  disks  for  later  reduction  and 
evaluation. 
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D.   SYSTEM  MODIFICATIONS 

For  some  of  the  early  data  taken  near  atmospheric  pres- 
sure, the  measured  steam-side  heat- transfer  coefficient, 
when  operating  with  a  closed  system,  disagreed  by  as  much  as 
±30  percent  between  different  data  runs.   This  was  later 
found  to  be  due  to  the  malfunctioning  of  the  relief  valve 
provided  on  the  apparatus  for  safety  reasons  (see  Figure 
2.1).   This  valve  had  a  designed  relief  gauge  pressure  of 
52  mmHg.   However,  occasionally  it  was  found  to  open  with  a 
system  gauge  pressure  of  only  about  7  mmHg.   Once  opened,  the 
valve  remained  open  even  if  the  system  operating  pressure 
was  10  mmHg  below  atmospheric,  thus  allowing  air  to  enter 
the  apparatus.   To  avoid  this  problem,  subsequent  runs  were 
made  without  the  relief  valve  and  a  new,  more  reliable, 
valve  was  ordered.   In  the  meantime,  extreme  care  had  to  be 
exercised  to  avoid  system  pressurization. 

Originally,  the  apparatus  was  provided  with  a  10.6  kW 
(3  ton)  air-conditioning  system  to  cool  the  cooling  water 
sump.   However,  as  also  mentioned  by  Poole  [Ref.  13],  this 
system  was  inadequate  to  keep  the  sump  temperature  from 
rising.   For  this  reason,  this  air-conditioning  system  was 
removed.   By  removing  this  unit,  the  cooling  water  friction 
losses  were  reduced  and  the  maximum  water  velocity  was  in- 
creased from  4.2  to  4.4  m/s .   In  the  previous  stage  of  the 
investigation,  fresh  tap  water  was  continuously  fed  to  a  sump 
while  an  equal  amount  of  water  was  being  drained,  maintaining 
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a  constant  sump  level.   Since  the  sump  was  small,  the  hot 
water  from  the  test  section  had  to  be  returned  to  the  sump. 
For  this  reason,  there  was  a  continuous  variation  of  the 
sump  water  temperature.   In  order  to  solve  this  problem,  the 
sump  was  replaced  by  an  aluminum  tank  three  times  its  size 
and  the  outlet  hot  water  was  drained  directly  after  its 
exit  from  the  test  tube. 

The  air  ejector  used  to  create  vacuum  in  the  test  sec- 
tion could  not  provide  a  successful  operation  owing  to  the 
cycling  of  the  feeding  air  pressure.   Occasionally  the  pres- 
sure in  the  test  section  could  become  smaller  than  the  pres- 
sure in  the  ejector's  air  line  resulting  in  a  back  flow 
which  could  contaminate  the  apparatus.   This  problem  was 
solved  by  the  use  of  a  vacuum  pump.   The  suction  line  for 
this  vacuum  pump  was  a  penetration  in  the  auxiliary  condenser. 
Unavoidably,  there  was  a  continuous  removal  of  steam  together 
with  the  non-condensing  gases.   To  avoid  steam  from  reaching 
the  vacuum  pump  (vacuum  pumps  are  designed  to  remove  only 
air) ,  the  mixture  was  passed  through  a  heat  exchanger  where 
the  steam  was  condensed.   This  heat  exchanger  was  installed 
as  shown  in  Figure  2.3.   Note  that  no  new  supply  of  cooling 
water  was  used  for  this  heat  exchanger;  instead,  the  tap  water 
supply  line  into  the  water  sump  was  used  for  this  purpose. 

Careful  and  persistent  investigations  into  the  occurrence 
of  partial  dropwise  condensation  on  the  copper  tube  revealed 
that  the  presence  of  nylon  components,  such  as  tube  bushings 
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(see  Figure  2.2)  and  numerous  ferrules  used  in  the  apparatus, 
were  a  source  of  system  contamination.   This  problem  worsened 
when  the  system  pressure  was  reduced  for  a  vacuum  run  after 
the  system  had  operated  at  near-atmospheric  pressures.   This 
was  thought  to  be  owing  to  the  fact  that  outgassing  increases 
with  increasing  temperatures  and  decreasing  pressures.   To 
solve  these  problems,  new  Teflon  tube  bushings  were  manufac- 
tured.  Also,  all  nylon  ferrules  were  replaced  with  Teflon 
ferrules.   Note  that  Teflon  has  an  outgassing  rate  about  one- 
hundred  times  smaller  than  nylon  [Ref.  14]. 

E.   VACUUM  INTEGRITY 

Vacuum  tightness  for  any  condensation  experiment, 
especially  under  low  pressure  conditions,  is  of  vital  impor- 
tance.  Any  data  taken  while  the  apparatus  is  inadequately 
air-tight,  will  lead  to  false  conclusions  because  of  the 
deleterious  effect  of  the  non-condensing  gases."  Many  efforts 
were  made  in  the  past  to  achieve  a  leak-free  system.   During 
this  thesis  effort,  there  were  leak  tests  on  a  routine  basis. 
A  total  of  eight  vacuum-tightness  tests  were  carried  out  over 
a  period  of  ten  months,  including  the  period  devoted  to  taking 
data.   The  leak  rate  was  always  less  than  that  represented 
by  a  pressure  rise  equal  to  2  mmHg  in  a  24-hour  period. 
These  tests  were  carried  out  at  a  pressure  of  about  85  mmHg 
absolute . 
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F.   TUBES  TESTED 

An  instrumented  tube  was  fabricated  from  a  thick-walled 
copper  tube  with  an  inner  diameter  of  12.7  mm  (1/2  in)  and 
an  outer  diameter  of  18  mm  (3/4  in).   The  condensing  length 
was  13.3  mm.   This  tube  was  provided  with  six  wall  thermo- 
couples, equally  spaced  around  the  tube  periphery  and 
located  radially  at  the  mid-point  of  the  tube  wall  thickness. 
Details  about  the  fabrication  and  instrumentation  of  this 
tube  are  given  by  Poole  [Ref.  13].   Figure  2.4  shows  a  photo- 
graph of  the  instrumented  tube.   This  tube  was  used  to 
determine  the  leading  constant  of  a  Sieder-Tate  type  equation 
for  the  inside  heat-transfer  coefficient. 

A  smooth  tube  of  copper  with  dimensions  similar  to  the 
instrumented  tube  was  also  fabricated.   The  purpose  of  the 
smooth  tube  was  to  determine  the  Sieder-Tate  leading  constant 
by  using  a  modified  Wilson  plot,  and  to  compare  the  data  for 
the  steam-side  coefficient  with  the  Nusselt  prediction. 

Twenty  one  finned  tubes  were  also  fabricated.   These 
tubes  had  the  same  inside  diameter  and  length  dimensions  as 
the  smooth  tube.   Further,  they  had  a  root  diameter  equal  to 
the  outside  diameter  of  the  smooth  tube.   The  fins  were 
rectangular  in  shape.   The  fin  dimensions  were  systematically 
changed  as  shown  in  Table  I.   Figure  2.5  shows  a  photograph 
of  the  group  of  tubes  with  a  fin  thickness  of  0.75  mm  and  a 
fin  height  of  1.0  mm. 
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TABLE  I 
Geometry  of  Finned  Tubes 


Fin 
Fin  Spacing  Thickness  Height 

(mm)  (mm)  (mm) 


0.5,  1.0,  1.5,  2.0,  4.0        1.0  1.0 

0.5,1.0,1.5,2.0,4.0        0.75  1.0 


0.5,1.0,1.5,2.0,4.0  0.5  1.0 

1.0,  1.5,  2.0,  4.0  1.0  2.0 

1.0  1.5  1.0 

1.25  1.25  1.0 


Th 

Fin 
ickness 
(mm) 

1 

.0 

0 

.75 

0, 

.5 

1, 

.0 

1, 

.5 

1, 

.25 

G.   INSERT 

Early  investigations  showed  that  with  no  internal  en- 
hancement, the  water-side  thermal  resistance  could  be  as  much 
as  90  percent  of  the  overall  resistance.   Under  such  condi- 
tions, even  a  small  discrepancy  on  the  water-side  resistance 
will  translate  into  a  substantially  greater  discrepancy  on 
the  condensing  heat-transfer  coefficient.   Therefore,  in 
order  to  improve  the  accuracy  of  the  condensing  heat-transfer 
coefficient,  a  decision  was  made  to  enhance  the  inside  heat- 
transfer  coefficient. 

For  this  purpose,  a  spiral  insert  was  fabricated.   The 
insert  consisted  of  a  6.4  mm  diameter  stainless-steel  rod 
together  with  a  copper  wire  that  was  wrapped  around  the  rod. 
The  diameter  and  the  pitch  of  the  wrapped  wire  were  3.2  mm 
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and  20  mm  respectively.   After  wrapping,  the  wire  was  properly 
soldered  on  the  rod  and  then  it  was  machined  down  to  give  a 
clearance  of  0.5  mm  between  the  outer  wire  diameter  and  the 
tube  inside  wall.   This  clearance  was  considered  important 
to  avoid  heat  conduction  from  the  tube  wall  to  the  insert. 
Figure  2.6  is  a  photograph  of  the  insert  with  one  of  the 
finned  tubes . 
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III.   SYSTEM  OPERATION  AND  PROBLEMS 

Each  tube  to  be  tested  was  cleaned  following  the  proce- 
dures listed  in  Appendix  A.   After  installing  a  tube,  the 
system  was  brought  to  operating  pressure  by  following  the 
procedures  listed  in  Appendix  B.   The  data  were  collected 
when  steady-state  conditions  were  achieved  which  were 
assumed  to  occur  when  the  temperature  of  the  steam  in  the 
test  section  and  the  inlet-outlet  water  temperatures  were 
stabilized.   Double  data  sets  were  taken  for  each  flow 
setting  starting  with  a  cooling  water  flow  rate  of  80% 
(4.4  m/s).   The  next  water  flowmeter  settings  were  70,  60, 
45,  35,  26,  and  20  percent  followed  by  55  and  80  percent. 
These  settings  were  selected  to  provide  nearly  equal-spaced 
heat  flux  data  points. 

Since  the  tubes  were  fabricated  from  copper,  a  highly 
conductive  material,  the  time  interval  required  to  achieve 
steady-state  conditions  was  not  more  than  two  minutes.   A 
whole  run  (18  data  sets)  was  completed  in  a  time  interval  of 
less  than  one  hour.   After  each  change  in  the  cooling  water 
flow  rate,  the  system  pressure  changed  slightly.   In  order 
to  keep  this  system  pressure  nearly  constant,  it  was  neces- 
sary to  adjust  the  flow  rate  through  the  auxiliary  condenser 

During  the  early  stages  of  this  experiment,  it  was 
observed  that  the  thermocouple  readings  were  disturbed  by  an 
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external  radio-frequency  (RF)  signal.   This  RF  signal  was 
generated  by  the  special  device  used  to  control  the  power 
input  to  the  boiler  heaters.   This  device  consisted  of  a 
silicon-controlled  rectifier  that  chopped  the  input  signal 
to  obtain  the  necessary  electrical  power,  while  the  chopped 
portion  was  emitted  as  an  RF  signal.   This  RF  signal  caused, 
for  example,  the  steam  temperature  to  fluctuate  within  ±4  K, 
while  the  thermopile  readings  fluctuated  within  ±0.5  K  (note 
that  the  thermopile  readings  measuring  temperature  rise  of 
cooling  water  can  be  as  low  as  0.5  K) .   To  alleviate  these 
problems,  all  thermocouples  were  shielded  against  the  RF 
signal  using  braided  copper  cables.   This  substantially 
reduced  the  temperature  fluctuations.   For  example,  after 
the  shielding  was  installed,  the  thermopile  readings  fluctu- 
ated less  than  ±0.04  K.   Further  efforts  revealed  that  the 
small  fluctuations  which  persisted  were  due  to  the  presence 
of  a  cable  that  carried  an  indication  of  the  voltage  to  the 
boiler  heaters.   This  voltage  is  of  "DC"  type  and  it  was  fed 
to  the  data  acquisition  system.   Unfortunately,  this  voltage 
was  not  truly  DC;  instead  it  consisted  of  the  RF  disturbances 
described  earlier.   When  this  cable  was  disconnected,  abso- 
lutely no  fluctuations  were  observed.   Thus,  a  switch  was 
added  on  this  cable  to  disconnect  the  source  of  the  RF 
signal  while  thermocouple  readings  were  being  recorded.   Soft- 
ware was  modified  to  gather  all  thermocouple  readings  while 
this  switch  was  turned  off.   When  the  computer  was  ready  to 
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'read'  the  voltage  to  the  boiler  heaters,  the  program  was 
coded  to  'pause'  and  then  remind  the  operator  to  turn  on  the 
switch.   Immediately  after  taking  this  reading,  the  computer 
reminds  the  operator  to  turn  off  this  switch  again. 

The  appearance  of  the  condensate  film  was  regularly 
observed  through  the  viewport.   In  the  event  the  condensate 
film  appeared  to  be  patchy,  the  data  run  was  discontinued 
and  the  data  already  collected  were  discarded.   During  some 
runs,  the  film  appeared  to  be  totally  filmwise,  but  the 
final  data  set  indicated  a  condensing  coefficient  up  to  10 
percent  greater  than  the  value  for  the  initial  data  point. 
Recall  that  the  initial  and  final  data  sets  were  for  the 
same  cooling  water  flow  rate.   It  is  worth  noticing  that  it 
was  easier  to  obtain  filmwise  condensation  under  vacuum 
conditions  than  at  near-atmospheric  pressures.   This  is 
probably  due  to  lower  vapor  velocities  associated  with 
atmospheric  runs  (2  m/s  for  vacuum  and  1  m/s  for  atmospheric 
conditions) .   All  the  data  presented  in  this  thesis  showed 
less  than  ±5  percent  disagreement  in  condensing  coefficient 
between  initial  and  final  data  sets. 

A.   THE  DROPWISE  CONDENSATION  PROBLEM 

The  main  objective  of  this  thesis  effort  was  to  take 
data  under  strict  filmwise  condensation  conditions  using 
copper  tubes.   Pure  copper  has  low  wettability  with  water  in 
the  presence  of  a  minute  contamination;  thus,  it  has  the 
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tendency  to  condense  steam  under  the  dropwise  mode.   Such  a 
minute  contamination  is  generally  difficult  to  avoid.   A 
slightly  hydrophobic  surface  condition  can  lead  to  very 
false  conclusions  since  the  dropwise  condensation  mode  of 
heat  transfer  is  much  more  effective  than  the  filmwise  mode. 

In  the  first  stage  of  this  thesis  effort,  the  tubes  were 
prepared  following  the  procedure  recommended  by  Poole 
[Ref.  13].   This  procedure  included  the  cleaning  of  tubes 
with  a  mixture  of  sodium  hydroxide  and  ethyl  alcohol  (equal 
amounts  by  mass)  at  a  temperature  of  about  80  °C.   This  proce- 
dure was  soon  proven  inadequate  as  the  copper  tube  showed 
partial  dropwise  condensation  within  20  minutes.   This  was 
verified  during  a  four-hour  endurance  test.   The  water  flow 
rate  was  held  constant  and  every  30  minutes  the  temperature 
rise  of  the  cooling  water  was  recorded.   This  rise  continu- 
ously increased  and  gradually  the  condensation  changed  to 
the  dropwise  mode. 

As  a  solution  to  this  problem,  it  was  decided  to  oxidize 
the  tube  surface  more  thoroughly  using  the  above-mentioned 
mixture  (50  percent)  aqueous  sodium  hydroxide  mixed  with  an 
equal  part  of  ethyl  alcohol) .   The  tube  was  painted  with  this 
mixture  and  was  then  heated  with  steam  for  about  an  hour 
with  no  water  flowing  inside  the  tube.   This  process  resulted 
in  the  formation  of  a  dark  oxide  layer  on  the  copper  surface, 
which  had  high  wetting  characteristics  and  an  insignificant 
thermal  resistance.   After  the  oxidation,  new  endurance  tests 
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were  made.   For  a  time  period  of  three  hours,  the  cooling 
water  temperature  rise  was  essentially  unchanged,  while  the 
condensate  film  appearance  was  excellent.   Steam-side  coeffi- 
cients obtained  with  and  without  oxidation  agreed  to  within 
±2  percent.   This  was  probably  due  to  the  avoidance  of  even 
minute  areas  undergoing  dropwise  condensation. 

The  most  significant  achievement  of  this  oxidation  was 
the  reduction  of  data  scattering  (less  than  ±1  percent) ,  and 
the  very  good  repeatability  of  similar  runs  (within  ±3  percent) 

B.   VAPOR  VELOCITY  LIMITATIONS 

The  major  assumption  in  the  Nusselt  theory  is  that  the 
vapor  velocity  does  not  induce  shear  forces  on  the  vapor- 
liquid  interface.   This  assumption  is  not  realistic  in  prac- 
tice.  In  order  to  regulate  the  pressure  in  the  test  section, 
it  is  necessary  to  have  some  steam  condensed  in  the  auxiliary 
condenser.   This  means  that  there  must  be  steam  flow  into 
the  auxiliary  condenser.   The  steam  flow  rate  (i.e.,  the 
vapor  velocity)  is  a  function  of  the  power  input  to  the 
heaters.   In  order  to  achieve  low  vapor  velocities,  the  power 
input  should  be  decreased.   At  low  pressures,  however,  it 
was  evident  that  the  fluctuation  of  system  pressure  increased 
with  a  decrease  in  boiler  power.   This  fluctuation  was  due 
to  the  intermittent  breakup  of  vapor  bubbles  on  the  water- 
vapor  interface  in  the  boiler.   Surprisingly,  this  small 
fluctuation  in  system  pressure,  and  hence  vapor  velocity, 
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reflected  in  the  cooling  water  outlet  temperature.   To 
minimize  this  fluctuation,  a  vapor  velocity  of  approximately 
2  m/s  past  the  tube  was  selected  for  vacuum  runs.   For  the 
atmospheric  runs  (maximum  power  input) ,  the  cooling  water 
outlet  temperature  fluctuation  was  less  than  ±0.005  K,  while 
for  the  vacuum  runs  (half  of  maximum  power  input) ,  the  fluc- 
tuation was  less  than  ±0.05  K. 

A  low  vapor  velocity  is  also  undesirable  from  the  film- 
wise  point  of  view.   When  the  vapor  flow  is  low,  then  the 
condensation  rate  decreases  and  the  film  can  very  easily 
break,  resulting  in  a  change  of  the  condensation  mode.   Once 
the  film  has  broken,  there  is  no  possibility  for  the  tube  to 
resume  the  filmwise  mode.   Also,  when  the  vapor  velocity  is 
high,  there  is  a  very  strong  wash-off  effect  on  the  tube, 
resulting  in  the  reduction  in  the  tendency  of  the  tube  to 
change  the  condensation  mode. 

It  was  decided  to  use  2  m/s  and  1  m/s  vapor  velocities 
for  the  vacuum  and  atmospheric  runs  respectively.   These 
velocities  are  considered  much  smaller  than  the  actual  ones 
in  condensers.   However,  due  to  the  limited  power  input,  the 
1.0  m/s  used  for  the  atmospheric  runs  was  indeed  the  maximum 
possible.   For  the  vacuum  runs,  the  2.0  m/s  vapor  velocity 
was  achieved  by  using  half  of  the  available  power.   Some 
attempts  were  made  to  increase  the  vapor  velocity  to  4.0 
m/s,  but  it  was  considered  unsafe  due  to  the  violent  and 
rapid  explosions  of  the  vapor  bubbles  in  the  boiler. 
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IV.   DATA  REDUCTION 

During  this  thesis  effort,  all  the  programs,  including 
property  functions  and  the  calibration  curves  for  the 
cooling  water  flowmeter  and  for  all  thermocouples,  were 
carefully  checked.   Only  a  very  few  minor  corrections  had  to 
be  made.   The  programs  have  the  following  noteworthy  features 

1.  In  order  to  obtain  a  meaningful  mixing-cup  temperature 
at  the  coolant  outlet,  a  mixing  chamber  was  provided  at  the 
outlet  as  shown  in  Figure  2.2.   As  also  discussed  by  Poole 

[Ref.  13],  the  considerable  pressure  drop  created  by  this 
mixing  chamber  introduced  a  temperature  rise  in  the  cooling 
water.   With  no  steam  condensation  taking  place,  the  coolant 
temperature  rise  was  measured  for  various  water  velocities. 
The  results  were  plotted  as  shown  in  Figure  4.1.   When  taking 
condensation  data,  the  coolant  temperature  rise  was  corrected 
by  subtracting  the  temperature  rise  due  to  this  frictional 
heating  at  the  mixing  chamber.   For  this  purpose,  the  least- 
squares  line  shown  in  Figure  4.1  was  used. 

2.  The  smooth  end  sections  of  the  tube,  which  were  neces- 
sary to  hold  the  tube  in  place  (see  Figure  2.2) ,  were  acting 
as  axial  fins.   Since  copper  has  a  very  large  thermal 
conductivity,  axial  heat  conduction  occurred.   This  heat 
conduction  was  incorporated  into  the  data  reduction  process 
by  computing  the  fin  efficiency  of  the  smooth  end  sections 
of  the  tube. 
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3.   The  tap  water  temperature  changes  considerably  depend- 
on  the  season:   it  is  about  15  °C  in  winter  and  about 
22  °C  in  summer.   Because  of  these  changes,  it  was  necessary 
to  incorporate  water  property  variations  in  calibrating  the 
flowmeter.   For  this  purpose,  the  following  function  was 
used,  as  done  so  in  [Ref .  15] ,  to  calculate  a  correction 
factor  for  viscosity  changes. 


Cf   =   1.0365-1.9644(10  3) Tin  + 5 . 25 ( 10  6)Tin2      (4.1) 


The  value  of  the  above  function  at  the  time  of  flowmeter 
calibration  (i.e.,  at  Tin  =  22.0  °C)  was  0.995.   Thus,  in 
order  to  obtain  the  actual  mass  flow  rate  for  an  inlet 
temperature  of  Tin,  the  following  equation  was  used: 


m  =      ni      -.      7r — n  „  -  (4.2) 

act  calc    0.995 


The  major  steps  for  the  data-reduction  programs  are 
listed  below. 

A.   PROGRAM  SIEDER 

This  program  calculates  the  constants  of  the  Sieder-Tate 
type  equation  from  the  instrumented-tube  data. 

1.  Read  the  cooling  water  inlet  and  outlet  temperatures 
and  compute  the  average  temperature. 

2.  Read  the  six  wall  thermocouples  and  compute  the 
average  wall  temperature. 
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3.  Convert  the  water  volume  flowmeter  setting  to  mass 
flow  rate  by  dividing  the  volumetric  water  flow  rate  with 
the  water  velocity. 

4.  From  continuity,  compute  the  cooling  water  velocity. 

5.  Correct  the  outlet  water  temperature  for  the  mixing 
chamber  effect  and  then  re-compute  the  average  temperature. 

6.  Compute  the  amount  of  heat  transferred  to  the  cooling 
water. 

7.  Compute  the  wall  thermal  resistance  and  then  find  the 
inner  wall  temperature  (from  step  2,  the  average  temperature 
at  the  radial  location  of  the  thermocouples,  is  known) . 

8.  Compute  the  log-mean- temperature  difference. 

9.  Compute  the  efficiency  of  the  left  and  right  tube 
extensions. 

10.  Compute  the  inside  heat- transfer  coefficient. 

11.  Compute  the  Nusselt  number. 

12.  Compute  the  leading  constant  of  the  Sieder-Tate 
coefficient  for  this  data  set. 

13.  Repeat  steps  1  to  12  for  all  the  data  sets  of  the  run. 

14.  Make  a  least-squares  straight  line  fitting  of  the 
form: 


n     n  °'8  t,  1/3   /C>0*14  +  B  (4.3) 

Nu   =   C  Re     Pr  '         — 

^w 


The  C  is  the  unknown  leading  constant,  while  B  is  an  addi- 
tional constant  found  to  give  an  improved  fitting. 
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B.   MODIFIED  WILSON  PLOT  PROGRAM 

This  program  calculates  the  leading  constant  of  the 
Sieder-Tate  equation  from  the  smooth-tube  data. 

1.  Assume  a  value  for  the  leading  constant  of  the  Sieder- 
Tate  equation. 

2.  Make  the  steps  1,  3,  4  and  5  of  previous  program 
"SIEDER. " 

3.  Compute  the  heat  flux  from  the  specific  heat  of  water, 
the  mass  flow  rate,  the  cooling  water  temperature  difference 
and  the  outside  tube  surface. 

4.  Compute  the  log-mean-temperature  difference  between 
the  steam  temperature  and  the  inlet  and  outlet  cooling  water 
temperatures . 

5.  Compute  the  overall  heat-transfer  coefficient. 

6.  Assume  an  outer  tube  surface  temperature  in  order  to 
determine  the  film  temperature  for  the  evaluation  of  conden- 
sate properties. 

7.  Calculate  the  outside  heat-transfer  coefficient  by 
using  the  Nusselt  equation. 

8.  Compute  the  outer  tube  surface  temperature. 

9.  Iterate  between  steps  7  and  8  until  two  successive 
outer  surface  wall  temperatures  are  within  ±1  percent. 

10.  Assume  a  viscosity  correction  factor  for  the 
Sieder-Tate  equation. 

11.  Compute  the  fin  efficiency  of  the  extending  parts  of 
the  tube. 
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12.  Compute  the  inside  heat- transfer  coefficient. 

13.  Compute  the  inner  surface  temperature  in  order  to  find 
the  film  temperature  on  the  coolant  side. 

14.  Compute  the  viscosity  correction  factor  and  iterate 
between  steps  11  through  14  until  two  successive  viscosity 
correction  factors  are  within  ±1  percent. 

15.  Compute  the  new  Sieder-Tate  coefficient  and  iterate 
between  steps  2  through  15  until  two  successive  Sieder-Tate 
coefficients  are  within  ±0.5  percent. 

16.  Compute  the  Sieder-Tate  coefficient  for  all  the  data 
sets  of  the  run. 

17.  Make  a  least  square  straight  line  fitting  to  establish 
the  final  leading  constant  that  satisfies  all  the  data  sets. 

C.   MAIN  DATA  REDUCTION  PROGRAM  (DRP4) 

This  program  was  used  for  computing  the  condensing  heat- 
transfer  coefficient  for  both  finned  and  smooth  tubes. 

1.  Make  steps  1  through  5  of  program  "SIEDER. " 

2.  Compute  the  log-mean  temperature  difference  between 
the  steam  temperature,  and  the  inlet  and  outlet  cooling 
water  temperatures. 

3.  Compute  the  heat  flux. 

4.  Compute  the  overall  heat-transfer  coefficient. 

5.  Compute  the  wall  conduction  resistance. 

6.  Compute  the  inside  heat- transfer  coefficient. 

7.  Compute  the  condensing  heat-transfer  coefficient  by 
subtracting  the  wall  and  the  inside  resistances  from  the 
overall  resistance. 
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V .   RESULTS  AND  DISCUSSION 

Numerous  data  runs  were  made  using  the  procedures 
described  in  Chapter  III.   Each  tube  was  tested  at  least 
three  times,  on  different  days,  in  order  to  obtain  repeata- 
ble  results.   Data  were  taken  for  all  the  tubes  at  both 
near-atmospheric  and  vacuum  (85  mmHg,  1.7  psia)  conditions. 
Complete  filmwise  condensation  was  maintained  for  the  data 
runs  and  the  mass  concentration  of  non-condensing  gases  was 
held  between  0.0  and  -1.5  percent  during  all  the  testings. 
The  negative  value  was  indicative  of  either  a  slight  super- 
heat or  an  inaccurate  manometer  reading.   It  was  estimated 
that  an  error  of  ±0.2  K  in  steam  temperature  or  an  error  of 
1.0  mmHg  in  system  pressure  translates  into  approximately 
±1.2  percent  error  in  the  non-condensing  gas  concentration. 
Thus,  the  non-condensing  gas  concentration  was  essentially 
zero  to  within  the  accuracy  of  the  measurements. 

A.   INSIDE  HEAT-TRANSFER  COEFFICIENT 
1.   Instrumented  Tube  Results 

The  instrumented  tube  was  used  to  determine  the 
inside  heat-transfer  coefficient  mainly  owning  to  the  use  of 
the  insert  to  enhance  the  heat- transfer  performance.   By 
the  use  of  the  instrumented  tube,  the  wall  temperature  was 
determined  after  averaging  the  six  thermocouple  readings. 
Therefore,  the  main  assumption  in  the  calculation  of  the 
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inside  heat-transfer  coefficient  is  that  the  wall  tempera- 
ture profile  is  uniform.   However,  Figures  5.1  and  5.2  show 
the  experimentally-measured  temperature  profiles  around  the 
tube,  for  the  maximum  and  minimum  cooling  water  flow  rate 
settings,  for  the  vacuum  and  near-atmospheric  runs,  respec- 
tively.  It  is  obvious  that  the  wall  temperature  profile  is 
not  uniform;  therefore,  this  direct  averaging  procedure  may 
introduce  an  uncertainty  the  magnitude  of  which  is  unknown. 
However,  owing  to  the  time  constraint,  further  evaluation 
of  this  uncertainty  was  considered  to  be  beyond  the  scope 
of  this  investigation.   These  figures  also  show  least-squares- 
fit  cosine  curves  of  the  form: 


AT  ,  6  /  C    1  X 

—  =   1  -  a  cos  —  (5.1) 

AT 


It  is  worth  noticing  that  the  temperature  profile  corres- 
ponding to  the  high  heat  flux  is  flatter  than  the  profile 
corresponding  to  the  low  heat  flux.   Also,  the  temperature 
profile  of  the  near-atmospheric  run  is  flatter  than  the 
temperature  profile  of  the  vacuum  run.   It  is  concluded  that 
as  the  heat  flux  increases,  hi  increases,  hence  ho/hi 
decreases  and  the  temperature  profile  becomes  more  uniform; 
therefore,  the  assumption  of  an  isothermal  wall  becomes  more 
reasonable.   This  is  one  more  reason  why  the  internal 
enhancement  (use  of  insert)  improves  the  accuracy  of  the 
data. 
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Tables  II  and  III  give  the  values  of  the  amplitude  of  the 
cosine  fit  for  all  data  points  taken,  under  both  pressure 
conditions . 

TABLE  II 
Instrumented  Tube  Results  for  the  Vaccum  Run 


ho/hi 


Heat  Flux 
(MW/m2) 

0. 

209 

0. 

207 

0. 

199 

0. 

191 

0. 

184 

0. 

173 

0. 

160 

ho 
(kW/m2.K) 

10. 

89 

11. 

28 

11. 

45 

11. 

61 

11. 

76 

12. 

05 

12. 

09 

0.301 

0.135 

0.348 

0.142 

0.430 

0.158 

0.507 

0.168 

0.567 

0.178 

0.692 

0.185 

0.832 

0.202 

TABLE  III 
Instrumented  Tube  Results  for  the  Atmospheric  Run 


ho 
(kW/m2.K)        ho/hi 


9.18 
9.19 
9.25 
9.38 
9.49 
9.75 
10.05 


Heat  Flux 
( MW/m2 ) 

0. 

488 

0. 

486 

0. 

479 

0. 

462 

0. 

447 

0. 

417 

0. 

384 

0.201 

0.115 

0.221 

0.119 

0.249 

0.133 

0.289 

0.137 

0.346 

0.149 

0.437 

0.164 

0.585 

0.179 
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Also  the  same  tables  contain  additional  information 
about  the  heat  flux,  the  outside  heat-transfer  coefficient 
and  the  ratio  of  outside  to  inside  heat-transfer  coefficients 
for  all  the  cooling  water  flow  rate  settings. 

By  using  the  average  reading  of  all  thermocouples 
for  the  wall  temperature,  the  complete  form  of  the  Sieder- 
Tate-type  equation  was  determined  to  be: 


0  8    1/1   ^c  0  14 
Nu   =   0.0635  Reu    Pr  7  J  (_Si)u'-Lil  +  26.4  (5.2) 

u 
w 


The  constant  B  =  26.4  is  an  additional  constant  found  to 
improve  the  fitting.   Note  that  the  leading  coefficient  is 
about  2.5  times  greater  than  the  well-known  Sieder-Tate 
coefficient  of  0.027  [Ref.  16]  for  plain  tubes.   This  signi- 
ficant difference  is  mainly  due  to  the  coiled  insert  used 
to  boost  the  inside  heat-transfer  coefficient.   A  discussion 
of  the  advantages  and  disadvantages  of  the  coiled  insert  was 
given  in  Chapter  II,  Section  G.   Figure  5.3  shows  the  data 
points  for  both  pressure  conditions  together  with  a  fit  of 
equation  (5.2)  . 

It  is  worth  mentioning  that  the  use  of  thermocouples 
to  determine  the  wall  temperature  may  introduce  yet  another 
uncertainty  (in  addition  to  the  incorrect  averaging  mentioned 
earlier)  in  the  measurement  of  each  local  wall  temperature. 
This  is  due  to  the  completely  different  thermal  conductivity 
of  the  wall  area  in  the  immediate  vicinity  of  the  thermocouple 
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hot  junction.   For  example,  the  thermocouple  wires  are 
covered  with  plastic;  therefore,  they  present  a  low  thermal 
conductivity.   Also,  in  order  to  embed  the  thermocouple, 
there  must  be  sufficient  clearance,  which  certainly  presents 
an  additional  thermal  resistance  to  the  heat  flow. 

The  advantages  of  the  use  of  the  instrumented  tube 
are: 

1.  It  does  not  depend  on  other  theoretical  formulations 
that  may  have  unavoidable  errors  owing  to  the  assumptions 
made  to  formulate  the  problem  (for  example,  the  Modified 
Wilson  Plot  technique  uses  the  Nusselt-type  equation  for  the 
outside  heat-transfer  coefficient) . 

2.  It  is  not  too  sensitive  to  vapor  shear  effects  since 
these  effects  are  lumped  into  the  thermocouple  readings 
(i.e.,  the  vapor  shear  reduces  the  thickness  of  the  conden- 
sate layer;  therefore,  the  thermocouples  would  give  higher 
temperatures) . 

A  computer  program,  named  "SIEDER,"  was  written  to 
reduce  the  data  taken  on  the  instrumented  tube.   For  com- 
pleteness, the  raw  data  taken  on  the  instrumented  tube  for 
both  pressure  conditions  are  given  in  Appendix  C,  Table  VI. 
Hereinafter,  the  data  reduction  method  that  uses  the  instru- 
mented tube  water-side  heat-transfer  coefficient  will  be 
referred  to  as  the  "Direct"  method. 
2 .   Uninstrumented  Tube  Results 

The  uninstrumented  smooth-tube  data  were  used  to 
determine  the  inside  heat-transfer  coefficient  by  the  use  of 
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the  modified  Wilson  plot  technique.   Hereinafter,  this  tech- 
nique will  be  referred  to  as  the  "Wilson"  method.   Also  the 
smooth  tube  results  concerning  the  outside  heat-transfer 
coefficient  were  used  as  a  basis  for  comparison  with  the 
f inned-tube  results . 

The  analysis  procedures  used  for  the  Wilson  method 
is  the  same  as  that  given  in  [Ref .  17] .   While  a  Sieder- 
Tate-type  equation  was  used  to  represent  the  inside  heat- 
transfer  coefficient,  a  Nusselt  type  expression  was  used  for 
the  outside  coefficient  as  given  by: 

k^pf(p     -p    )h        1/3 

ho    =    a[ i^iyi ]  <5-3> 

where  a  is  a  disposable  constant  that  must  be  found  itera- 
tively.   Note  that  with  zero  vapor  velocity,  the  Nusselt 
prediction  results  in  a  value  of  0.655  for  a.   However, 
since  it  is  impossible  to  obtain  zero  vapor  velocity,  a 
higher  value  is  generally  obtained.   During  this  work,  an 
approximate  value  of  0.7  2  was  obtained  for  a. 

A  computer  program,  named  "WILSON,"  was  written  to 
run  the  Wilson  method  on  the  smooth-tube  data.   A  total  of 
four  runs  (two  each  under  vacuum  and  at  near-atmospheric 
pressures)  were  made.   These  runs  resulted  in  values  of  the 
Sieder-Tate  coefficient  of  0.0716,  0.0714,  0.0700  and 
0.070  8,  with  an  average  value  of  0.071.   This  value  is  about 
10  percent  higher  than  the  value  found  with  the  instrumented 
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tube;  i.e.,  the  Wilson  method  results  in  a  10  percent 
smaller  water-side  resistance  (less  conservative)  than  the 
resistance  predicted  by  the  use  of  the  instrumented  tube. 
Again  for  completeness,  the  raw  data  taken  on  the  uninstru- 
mented  tube  under  both  pressure  conditions  are  given  in 
Appendix  C,  Table  V. 

3 .   Summary  of  the  Inside  Heat-Transfer  Coefficient 

In  order  to  compare  the  two  above-mentioned  methods 
(Direct  and  Wilson) ,  the  data  taken  on  a  smooth  tube  under 
vacuum  conditions  were  reduced  to  yield  the  outside  heat- 
transfer  coefficients  using  the  water-side  coefficients 
derived  by  each  method.   For  this  purpose,  these  data  were 
plotted  in  Figure  5.4  together  with  the  original  Nusselt 
equation  for  a  single  horizontal  tube,  re-written  in  terms 
of  the  indicated  non-dimensionalized  terms: 


Nu  ,  /c- 

=   0.725  F  '  (5.4) 


Re0'5 


Also   shown    in   this    figure   is    the   correlation   of    Fujii,    et 
al    [Ref .    18] : 


Nus  1/5 

=      0.96    F   /D  (5.5) 


Re0-5 


The  Fujii  et  al.  correlation  is  based  on  numerous  data 
taken  on  tube  bundles  in  both  in-line  and  staggered  arrange- 
ments.  It  can  be  seen  that  the  coefficient  on  the  right-hand 
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side  of  equation  (5.5)  is  larger  than  the  value  in  equation 
(5.4)  owing  to  the  presence  of  vapor  shear  (recall  that  the 
Nusselt  equation  is  derived  for  zero-vapor-shear  conditions) . 
For  completeness,  the  error  bands  of  the  data  used  for  this 
correlation  are  also  shown. 

Note  in  Figure  5.4  that  the  data  of  this  investiga- 
tion on  the  smooth  tube  are  slightly  higher  than  those 
predicted  by  equation  (5.5)  .   Also,  it  is  worth  noting  that 
the  use  of  the  Wilson  method  results  in  data  that  are  in 
better  agreement  with  the  Fujii  correlation  than  the  data 
given  by  the  use  of  the  Direct  method.   Thus,  with  an  appro- 
priate function  for  the  outside  heat-transfer  coefficient 
and  accurate  smooth-tube  data,  the  modified  Wilson  plot 
technique  appears  to  offer  a  satisfactory  and  simple  method 
for  obtaining  the  inside  heat-transfer  correlation.   However, 
since  the  error  band  of  the  Fujii  correlation  is  much  wider 
than  the  disagreement  between  this  correlation  and  the 
present  data,  it  may  be  unwise  to  draw  any  conclusions  as 
to  the  accuracy  of  these  two  methods  (Direct  or  Wilson) . 
As  a  result,  it  is  not  possible  to  say  which  method  of  measur- 
ing the  inside  heat-transfer  coefficient  is  more  accurate. 
However,  except  where  noted  below,  all  the  outside  heat- 
transfer  coefficients  obtained  in  this  thesis  were  calcu- 
lated using  the  Direct  method. 
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B.   PERFORMANCE  OF  FINNED  TUBES 

Data  were  taken  for  the  finned  tubes  given  by  Table  I, 
under  vacuum  (85  mmHg)  and  at  near-atmospheric  pressures 
with  vapor  velocities  of  2  m/s  and  1  m/s  respectively.   The 
computed  steam-side  coefficient  was  based  on  the  smooth  tube 
outside  area  (i.e.,  the  area  if  the  fins  were  machined  off). 

The  data  reductions  were  performed  using  the  program 
named  "DRP4."   This  is  a  revised  version  of  the  program 
(DRP2)  used  by  Poole  [Ref.  13].   The  raw  data  taken  on  finned 
tubes  are  presented  in  Appendix  D. 

1 .   Condensation  on  Fins  in  General 

Figure  5.5  shows  the  variation  of  heat  flux  with  the 
steam-side  temperature  drop  for  a  set  of  five  finned  tubes 
under  vacuum  conditions.   Both  the  thickness  and  height  of 
the  fins  were  1.0  mm,  while  the  spacings  were  0.5,  1.0,  1.5, 
2.0  and  4.0  mm.   Notice  that  the  slope  of  the  curves  in- 
creases with  fin  spacing  up  to  the  point  where  the  optimum 
spacing  (1.5  mm)  is  reached  and  then  the  slope  decreases 
with  further  increase  in  fin  spacing.   This  means  that  the 
enhancement  of  the  condensation  process  is  highly  dependent 
on  the  fin  spacing.   The  data  were  fitted  with  least-squares 
curves  using  the  following  relation: 


q   =   a  ATb  ,  (5.6) 


and  AT  was  computed  by  using  the  following  relation: 
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q   =   h   AT  (5.7) 


where  q  was  measured  experimentally,  while  h0  was  obtained 
by  subtracting  inside  and  wall  resistances  from  the  measured 
overall  resistance.   From  equation  (5.6),  it  is  possible  to 
calculate  the  heat  flux  for  any  temperature  difference  be- 
tween the  steam  and  the  tube  wall  when  the  constants  a  and 
b  are  known . 

Table  IV  gives  the  values  of  a  and  b  for  all  the 
tubes  tested  for  both  pressure  conditions.   Notice  that  the 
value  of  the  constant  a  has  almost  the  same  range  for  all 
the  finned  tested  tubes  for  both  pressure  conditions.   The 
values  of  both  constants  a  and  b  increase  as  the  fin  spacing 
increases  up  to  the  spacing  of  1.5  mm  (it  will  be  shown 
later  that  this  is  the  optimum  spacing)  and  then  they 
decrease.   The  range  of  the  exponent  b  is  different  between 
the  vacuum  and  near-atmospheric  runs;  i.e.,  the  range  for 
the  vacuum  runs  is  0.76  to  0.89,  while  for  the  near-atmospheric 
runs  it  is  0. 87  to  0.94. 

Figure  5.6  shows  a  relative  plot  of  the  outside 
condensing  coefficient  for  the  same  set  of  tubes  just 
described.   For  comparison  purposes,  the  smooth-tube  data 
and  the  Nusselt  prediction  are  plotted  on  the  same  figure. 
All  tubes  tested  show  good  enhancement  over  the  smooth  tube. 
It  is  worth  noticing  that  even  the  smooth  tube  results  are 
higher  than  the  Nusselt  prediction;  this  is  mainly 
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attributable  to  Nusselt's  conservative  assumption  which 
neglects  vapor  shear. 

Combining  equations  (5.6)  and  (5.7) ,  h   can  be 
written  as: 


h   =   a1/b  g<»-D/b  (5.8) 

o  ^ 


It  appears  that  equation  (5.8)  results  in  curves  almost 
parallel  to  the  Nusselt  prediction,  indicating  that  the 
effect  of  heat  flux  is  the  same  on  the  finned  surface  as  on 
the  smooth  tube. 

On  the  above  figure,  the  dotted  lines  represent  the 
fittings  of  the  data  when  the  data  were  processed  by  the 
Wilson  method.   Notice  that  the  condensing  coefficients 
given  by  the  Wilson  method  are  lower  by  about  10  percent 
than  those  given  by  the  Direct  method  (recall  that  the 
water-side  resistance  given  by  the  Direct  method  is  higher 
than  that  given  by  the  Wilson  method) . 

It  is  also  shown  in  the  above  figure  that  as  the 
heat  flux  decreases,  the  vapor-side  heat-transfer  coeffi- 
cient increases.   Note  that  there  is  about  a  10  percent 
difference  between  the  high  and  the  low  heat  flux  results. 
The  reason  for  this  increase  is  that,  as  the  flux  decreases, 
the  condensate  thickness  decreases,  resulting  in  less 
conduction  resistance  and  improved  heat-transfer  coefficient. 


74 


2 .   Significance  of  Fin  Spacing 

This  section  presents  results  showing  the  variation 
of  the  condensing  coefficient  with  heat  flux,  as  well  as  the 
variation  of  heat  flux  with  steam-side  temperature  drop  for 
both  pressure  conditions,  for  the  four  sets  of  finned  tubes. 
Figures  5.10  to  5.25  show  that  for  all  the  fin  thicknesses 
(i.e.,  1.0,  0.75  and  0.5  mm)  and  for  each  of  the  fin  heights 
(i.e.,  1.0  and  2.0  mm),  the  best  performance  was  obtained 
with  the  tubes  of  1.5  mm  fin  spacing.   This  was  true  for 
both  pressure  conditions.   The  second  best  performance  was 
given  by  the  tubes  with  spacings  of  either  1.0  mm  or  2.0  mm. 
The  third  best  performance  was  given  by  the  tubes  with  fin 
spacing  of  4.0  mm  and  the  tubes  with  a  fin  spacing  of  0.5  mm 
showed  the  poorest  performance.   The  performances  of  the 
last  two  tubes  however  were  interchanged  as  the  pressure 
conditions  were  changed  from  vacuum  to  near-atmospheric .   An 
effort  to  explain  this  behavior  will  be  made  in  the  next 
sections,  where  the  influence  of  the  fin  thickness  and 
height  on  the  heat-transfer  performance  is  discussed. 

As  shown  in  Appendix  E  (Error  Analysis) ,  the  error 
in  the  calculation  of  the  condensing  coefficient  increases 
as  the  heat  flux  decreases.   For  this  reason,  the  comparison 
of  finned  tubes  should  be  done  in  the  regime  where  the  error 
is  the  minimum  possible.   This  means  that  the  data  points 
corresponding  to  the  high-heat-flux  settings  should  be 
chosen.   The  comparison  of  the  finned  tubes  is  made  through 
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the  enhancement  ratio,  which  is  defined  to  be  the  ratio  of 
the  outside  heat-transfer  coefficient  of  a  finned  tube  to 
the  smooth-tube  value  (both  evaluated  on  the  same  heat 
flux) .   However,  the  capability  of  the  smooth  tube  to  remove 
heat  is  limited  (see  Figure  5.6  showing  that  the  maximum 
heat  flux  on  the  smooth  tube  is  only  about  30  percent  of  the 
maximum  heat  flux  on  the  best  finned  tube) .   Thus,  it  was 
necessary  to  extrapolate  the  smooth-tube  data  to  higher  heat 
fluxes  in  order  to  make  the  comparison  of  the  finned  tubes 

on  an  equal  heat  flux  basis. 

2 
Choosing  0.35  and  1.0  MW/m   as  reference  heat  fluxes 

for  the  vacuum  and  near-atmospheric  runs,  respectively, 
cross  plots  of  the  enhancement  ratio  versus  the  fin  spacing 
are  given  in  Figure  5.7  for  the  set  of  tubes  with  1.0  mm  fin 
thickness  and  height.   It  is  shown  that  the  optimum  fin 
spacing  found  is  1.5  mm.   The  0 . 5-mm  fin  spacing  which  pro- 
vides the  most  area  increment  has  one  of  the  worst  performances. 
This  poor  performance  can  be  easily  explained  by  the  fact 
that  the  entire  interfin  area  was  covered  by  the  condensate. 
Visual  observations  showed  that  the  condensate  retention 
angle  was  180  degrees,  which  is  also  very  well  predicted  by 
equation  (1.12).   The  thermal  resistance  induced  by  this 
thick  condensate  layer  deteriorates  any  benefit  from  in- 
creased area;  therefore,  the  heat-transfer  performance  is 
highly  reduced. 
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As  the  fin  spacing  increased  from  0.5  to  1.0  mm, 
the  retention  angle  decreased  and  a  miniscus  in  the  interfin 
area  at  the  top  of  the  tube  appeared  (see  Figure  5.8(a)) . 
For  the  1.0  mm  fin  spacing,  the  retention  angle  was  about 
110  degrees.   The  change  of  the  flooding  characteristics  is 
the  reason  for  the  sudden  increase  of  the  enhancement  ratio 
as  the  fin  spacing  changed  from  0.5  to  1  mm.   Visual  obser- 
vations showed  that  the  optimum  fin  spacing  (1.5  mm)  resulted 
in  a  condensate  film  profile  (in  the  interfin  area  at  the 
top  part  of  the  tube)  similar  to  that  shown  in  Figure  5.8(b) . 
Notice  that  the  bottom  of  the  miniscus  is  just  touching  the 
tube  surface.   When  the  fin  spacing  is  further  increased, 
the  miniscus  radius  becomes  larger  and  the  retention  angle 
reduces  even  more.   At  the  2.0  mm  fin  spacing,  the  miniscus 
is  no  longer  an  arc  of  a  circle  (see  Figure  5.8(c)). 

As  mentioned  in  the  beginning  of  this  chapter,  the 
condensing-coef f icient  calculation  is  based  on  the  smooth 
tube  external  diameter;  therefore,  the  effect  of  area  addi- 
tion is  lumped  together  with  all  the  other  factors  that  may 
influence  the  heat- transfer  process.   In  order  to  give  some 
insight  to  the  effect  of  all  the  other  factors,  the  enhance- 
ment ratio  was  divided  by  the  area  gain  (i.e.,  the  ratio  of 
the  condensing  area  of  finned  tube  to  the  area  of  smooth 
tube)  of  each  finned  tube  in  order  to  obtain  the  normalized 
enhancement  ratio.   Figure  5.9  shows  this  normalized  enhance- 
ment ratio.   It  is  clear  that  the  enhancement  is  due  to 
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s  =  1 . 0  mm 
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(a) 


(b) 


s  =  2 . 0  mm 


(c) 


Figure  5.8.   Miniscus  Profile  for  Various  Fin  Spacings 

(a)  s  =  1.0  mm,  (b)  s  =  1.5  mm,  (c)  s  =  2.0  mm 


79 


CD 


ID 


(T) 


(\l 


ID 
CO 


ro 


ID 


ID 


oj 


ID 


C 

0   1 
fd 

en    • 


CD 


CO 

0  TS 

to  g 

M    <d 

0 

>  o 

^ 

0  ^ 

£ 

•H 

\ 

-P    II 

CJi 

(13 

C 

Pi  -p 

•H 

O 

4->  m 

ret 

G  0 

U. 

cd 

cn 

e  co 

cd  c 

c 

U  -H 

•H 

C  fa 

fa 

rd 

jC  -C 

C    4-> 

fa  -H 

5 

T3 

CD     CO 

N    CD 

■H  X! 

rH     d 

(d  E-" 

g 

u  u 

0   0 

2  M-l 

in 

cd 
u 

0 

& 

■H 

fa 


ot^ph   q.u8uiaouequa   pszT-^uiiON 


80 


more  than  just  an  area  gain.   This  verifies  the  measure- 
ments of  Yau  et  al.  [Ref .  3] ,  and  an  explanation  for  this 
behavior  was  presented  in  Chapter  I.   Also,  Figure  5.9  shows 
that  the  normalized  enhancement  ratio  increases  as  the 
spacing  increases  and  that  the  rate  of  this  effect  gets 
smaller  for  fin  spacings  larger  than  1.5  mm.   This  reduction 
is  attributable  to  the  fact  that  the  pulling  action  (thinning 
effect)  diminishes  as  the  fins  get  further  apart.   Equation 
(5.9)  gives  that  the  pressure  gradient  between  points  C  and 
B  (see  Figure  1.1)  is  inversely  proportional  to  the  increase 
of  fin  spacing  (s). 


dPx      „      APCB 


dx  1         s/2 
'  CB 


Therefore,  the  thinning  effect  decreases  as  the  fin  spacing 
increases.   It.  can  be  concluded  that  any  assumption  that 
neglects  the  surface-tension  effect  may  not  be  realistic 
when  modelling  the  steam  condensation  process  on  fins . 

Careful  observation  of  Figure  5.7  shows  that  the 
performance  of  the  tubes  of  0.  5  and  4.0-mm  fin  spacing 
interchanges  between  the  vacuum  and  atmospheric  runs, 
respectively.   The  heat  that  was  removed  during  the  vacuum 
runs  was  about  three  times  smaller  than  the  heat  removed 
during  the  near-atmospheric  runs;  thus,  the  rate  of  conden- 
sate generation  for  the  vacuum  runs  was  three  times  smaller 
than  for  atmospheric  runs.   For  all  the  vacuum  runs,  no 
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matter  which  tube  was  used,  the  condensate  flow  pattern  was 
invisible.   During  the  near-atmospheric  runs,  the  flow 
pattern  was  also  invisible  except  for  the  tubes  of  0.5  mm 
fin  spacing.   In  these  cases  (0.5  mm  fin  spacing),  the 
draining  rate  was  very  fast  and  ripples  in  the  interfin 
space  were  clearly  visible.   If  this  is  the  case,  then  the 
convective  action  of  the  condensate  layer  should  enhance 
the  condensation  process.   After  this  observation,  it  became 
clear  why  the  performance  of  the  tubes  with  0.5  mm  spacing 
became  better  than  the  performance  of  the  tubes  with  the 
4.0  mm  spacing,  for  the  near-atmospheric  runs. 

As  mentioned  earlier,  Honda  and  Nozu  [Ref.  9] 
present  a  detailed  comparison  of  various  theoretical  models 
to  predict  the  condensing  coefficient  on  finned  tubes.   All 
the  models  have  the  assumption  that  the  inertia  term  in  the 
momentum  equation  and  the  convective  term  in  the  energy 
equation  can  be  neglected  in  comparison  with  the  other 
terms.   The  author  feels  that  the  above  assumption  is 
realistic  for  all  the  fin  spacings  except  those  that  corres- 
pond to  completely  flooded  tubes.   To  further  support  this 
argument,  the  reader  is  referred  to  [Ref.  9]  where  comparison 
of  the  four  models  is  presented  with  various  experimental 
data,  among  which  are  the  steam  data  given  by  Rose  [Ref.  3]. 
Honda  and  Nozu  show  that  the  maximum  disagreement  between 
the  data  of  Rose  and  predictions  (all  four  models)  is  for 
the  fin  spacing  of  0.5  mm  (completely  flooded  tube) . 
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3 .   Significance  of  Fin  Thickness 

Thcree  sets  of  finned  tubes  were  manufactured  with 
fin  thicknesses  of  0.5,  0.75  and  1.0  mm,  while  the  fin 
height  was  held  constant  at  1.0  mm.   Each  set  consisted  of 
five  tubes  with  fin  spacings  of  0.5,  1.0,  1.5,  2.0  and  4.0 
mm.   The  relative  plots  of  the  outside  heat- transfer  coeffi- 
cient  for  each  group  of  tubes,  tested  under  vacuum  and 
atmospheric  conditions,  were  already  shown  in  Figures  5.10 
to  5.21.   Notice  that  the  shape  of  the  curves  to  fit  the 
data  points  is  very  nearly  independent  of  the  fin  thickness. 
Figure  5.26  is  a  cross  plot  showing  the  enhancement  ratio  as 
a  function  of  the  fin  spacing  for  both  pressure  conditions, 
with  the  fin  thickness  as  a  parameter.   Notice  that  for 
almost  all  the  fin  spacings,  the  0 . 5-mm  fin  thickness  gave 
the  smallest  gain  even  though  it  provides  the  maximum 
condensing  area.   An  effort  to  explain  this  phenomenon  will 
be  made  at  the  end  of  this  section.   Between  the  0.75  and 
1.0-mm  fin  thicknesses,  it  is  not  very  clear  which  of  the 
two  is  superior.   For  the  vacuum  runs  with  fin  spacing 
around  the  optimum  value  (i.e.,  0.5  mm  <  s  <  1.5  mm)  it 
appears  that  the  1.0-mm  thickness  is  a  stronger  candidate 
for  the  optimum  fin  thickness.   On  the  other  hand,  for 
atmospheric  runs  with  most  fin  spacings,  the  0.75-mm  fin 
thickness  appears  to  be  the  best  candidate. 

In  order  to  further  study  the  dependence  of  fin 
thickness  on  the  enhancement  ratio,  data  shown  in  Figure 
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5.26  were  cross  plotted  as  shown  in  Figure  5.27.   As  mentioned 
earlier,  initially  there  were  only  three  tubes  with  fin 
thickness  of  0.5,  0.75  and  1.0  mm.   For  vacuum  conditions, 
the  inadequacy  of  this  fin  thickness  range  can  be  seen 
since  the  optimum  thickness  has  not  been  bracketed.   There- 
fore, an  additional  tube  with  a  fin  thickness  of  1.5  mm, 
a  fin  height  of  1.0  mm  and  a  fin  spacing  of  1.0  mm,  was 
manufactured.   The  enhancement  ratio  found  on  this  tube  is 
also  plotted  in  Figure  5.27.   It  can  be  seen  that  the  tube 
with  a  1.0-mm  fin  thickness  is  still  the  best  candidate  for 
the  vacuum  runs.   This  figure  also  shows  possible  error 
bands  for  the  computed  enhancement  ratios.   Since  the 
possible  error  bands  are  in  excess  of  the  enhancement  ratios 
at  these  two  fin  thicknesses,  no  attempts  were  made  to 
explain  why  the  optimum  fin  thickness  is  affected  by  the 
pressure  condition.   Therefore,  it  can  be  concluded  that  the 
optimum  fin  thickness  is  around  0.75  to  1.0  mm  for  both 
pressure  conditions. 

As  the  fin  thickness  increases  with  fin  spacing  held 
constant,  the  fin  density  (fins/mm)  decreases.   In  the 
limit,  the  tube  will  become  a  smooth  tube  with  an  outer 
diameter  qual  to  the  fin  diameter.   Recall  that  in  the 
Nusselt  equation  (5.3),  the  tube  diameter  is  in  the  denominator; 

therefore,  the  condensing  coefficient  is  inversely  propor- 

1/4 
tional  to  the  D     for  constant  wall  temperature,  or,  at 

o 

a  constant  wall  heat  flux,  the  condensing  coefficient  is 

1/3 
inversely  proportional  to  D    .   Thus, 
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D   1/3 

h:     =   h  (=r^)  (5.10) 

t+°°       O  D  ' 


where  h*    is  based  on  the  outer  diameter  of  the  tube, 
including  the  fin  height.   Now,  if  this  condensing  coeffi- 
cient is  based  on  root  diameter  (as  it  is  done  throughout 
this  thesis) ,  then 


h,.   D   =  h'   D^  (5.11) 


The  combination  of  equations  (5.10)  and  (5.11)  results  in 


Df  2/3 

h.     =   h  (=i)  (5.12) 

t->°°       o  D 
o 


This  means  that  there  will  be  very  little  enhancement  by 

increasing  the  fin  thickness  to  the  limit.   For  example, 

with  Dr  =  21.0  mm  and  D   =19.0  mm,  the  ratio  Iv   /h 
f  o  t-*00'  o 

would  equal  1.07. 

When  the  fin  thickness  is  reduced,  there  is  a  finite 
number  of  fins  that  can  be  added  in  a  given  tube  length 
(i.e.,  it  is  not  physically  possible  to  take  a  smooth  tube 
when  fins  of  infinitesimal  thickness  are  present) .   Since 
there  is  a  significant  increase  in  heat-transfer  area  as  the 
fins  become  thinner,  an  enhancement  for  the  condensing 
process  may  be  expected  on  the  basis  of  area  increment. 
Unfortunately,  this  was  not  the  case  as  discussed  below. 
Equation  (1.12)  predicts  the  retention  angle  and  it  is  shown 
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that  the  condensate  retention  is  independent  of  fin  thick- 
ness.  This  was  visually  verified  during  this  investigation. 
Comparing  two  tubes  which  have  the  same  length  and  fin 
spacing,  the  one  that  has  thinner  fins  would  hold  more 
condensate  when  the  retention  angle  is  the  same.   This  means 
that  part  of  the  metal  of  a  thick  fin  is  replaced  by  conden- 
sate, or  in  other  words,  the  thin-finned  tube  presents  a 
smaller  effective  thermal  conductivity  than  the  thick-finned 
tube  resulting  in  lower  performance.   Also,  as  the  fin 
thickness  decreases,  the  fin  efficiency  may  play  a  role  for 
tubes  made  of  low- thermal-conductivity  metal.   However,  for 
copper  tubes  with  fins  of  0.5  mm  thickness  and  1.0  mm 

height,  the  fin  efficiency  was  estimated  as  99  percent, 

2 
when  the  outside  heat-transfer  coefficient  was  40  kW/m  •  K. 

Since  both  fin  thickness  extremes  show  that  the 
heat-transfer  process  is  not  enhanced,  it  appears  that  the 
optimum  thickness  found  (around  0.75  to  1.0  mm)  is  reasonable 

Careful  observation  of  Figure  5.26  reveals  that  the 
tubes  with  0 . 5-mm  fin  spacing  do  not  follow  the  trend  of  the 
tubes  with  larger  fin  spacings.   In  this  case,  the  maximum 
performance  was  given  by  the  tube  with  the  0.75-mm  thick- 
ness, then  by  the  tube  with  0 . 5-mm  thickness  and  finally  by 
the  tube  with  1.0-mm  fin  thickness.   This  sequence  holds  for 
both  pressure  conditions.   The  rate  of  condensate  flow 
around  the  tube  was  observed  to  be  dependent  on  the  fin 
thickness.   The  fastest  flow  pattern  was  observed  for  the 
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tube  with  0.5  mm  fin  thickness,  then  for  the  tube  with  0.75 
mm  and  finally  for  the  tube  with  1.0  mm  fin  thickness.   This 
surprising  observation  can  be  explained  if  the  fin  tip 
effect  is  introduced;  i.e.,  since  the  fin  tip  is  rectangu- 
lar, there  is  a  possibility  that  some  of  the  condensate 
will  fall  into  the  interfin  spacing  and  some  will  follow  the 
circumferential  path.   If  this  condensate  flow  is  realistic, 
it  can  be  concluded  that  the  thicker  the  fin,  the  more 
condensate  will  flow  circumf erentially  along  the  fin  tip. 
To  further  support  this  conclusion,  numerous  runs 
were  made  at  a  different  fin  spacing  (1.0  mm)  where  the 
miniscus  in  the  interfin  area  at  the  top  part  of  the  tube 
is  visible  (recall  that  the  tubes  with  0 . 5-mm  fin  spacing 
are  completely  flooded) .   If  the  miniscus  radius  is  different 
as  the  fin  thickness  changes,  it  can  be  concluded  that  the 
condensate  flow  pattern  is  influenced  by  the  fin  thickness. 
The  tube  with  0.5  mm  fin  thickness  had  a  smaller  miniscus 
radius  than  the  tube  with  1.0  and  1.5  mm  fin  thickness. 
This  means  that  the  previous  conclusion  may  be  correct; 
therefore,  when  the  fin  tip  is  wide  and  of  course  rectangu- 
lar, there  will  be  circumferential  drainage  along  the  fin 
tip.   On  the  other  hand,  for  tubes  with  narrow  fin  tips, 
most  of  the  condensate  will  fall  into  the  interfin  spacing, 
resulting  in  a  very  fast  flow  pattern  and  the  creation  of 
ripples.   If  this  is  the  case,  the  thin  fins  should  show 
an  increased  performance  due  to  the  convective  action  of 
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the  condensate.   However,  the  final  performance  is  limited 
by  the  excessive  condensate  amount  that  is  attached  continu- 
ously on  the  tube.   These  may  be  the  reasons  why  the  tube 
with  0.5  mm  fin  thickness  has  better  performance  than  the 
tube  with  1.0  mm  fin  thickness  when  the  fin  spacing  is 
0.5  mm.   Between  the  tubes  with  0.5  and  0.75  mm  fin  thick- 
nesses, there  is  no  clear  visible  difference  in  the  flow 
pattern;  therefore,  the  last  tube  is  superior  because  it 
presents  less  conduction  resistance  (it  holds  less  conden- 
sate) ,  and  it  certainly  has  a  higher  fin  efficiency. 

In  general,  it  can  be  concluded  that  for  fins  of  1.0_ 
mm  height,  the  enhancement  ratio  was  weakly  dependent  on  fin 
thickness .   The  most  significant  reason  is  that  the  fin 
thickness  has  no  influence  on  the  retention  problem  as  the 
fin  spacing  does . 

In  order  to  check  the  validity  of  some  of  the 
previous  conclusions,  data  were  taken  for  one  more  tube. 
This  tube  had  a  fin  spacing,  thickness  and  height  of  1.25 
mm,  1.25  mm  and  1  mm,  respectively,  and  it  was  made  from  the 
tube  with  1.5  mm  fin  thickness  whose  fins  were  machined  down 
by  0.25  mm.   A  tube  with  1.25-mm  fin  spacing  was  expected 
to  have  a  better  heat-transfer  performance  than  the  tube  of 
1.0  mm  fin  spacing  since  the  condensing  heat-transfer  coeffi- 
cient increases  up  to  a  fin  spacing  of  around  1.5  mm. 
However,  the  performance  should  be  limited  owing  to  the 
fact  that  the  fin  thickness  is  outside  the  optimum  range 
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(0.75  to  1.0  mm).   Figures  5.28  and  5.29  show  the  relative 
plots  of  the  condensing  coefficient  with  heat  flux  for  both 
pressure  conditions,  for  three  tubes  whose  fin  geometry  is 
shown  in  these  figures.   It  is  shown  that  the  performance 
of  the  new  tube  is  almost  the  same  as  with  the  tube  with 
1.0  mm  fin  spacing  and  thickness.   This  means  that  what  was 
gained  by  the  spacing  increment  (from  1.0  to  1.25  mm)  was 
lost  by  the  fin  thickness  increase  (from  1.0  to  1.25  mm); 
therefore  the  previous  conclusions  about  the  fin  spacing  and 
thickness  effects  seem  to  be  reasonable. 
4 .   Significance  of  Fin  Height 

As  shown  in  Table  I,  two  fin  heights  (1.0  and  2.0 
mm)  were  examined  during  this  thesis  effort.   The  geometry 
of  tubes  tested  is  given  in  Table  I.   The  heat-transfer 
performance  of  these  tubes  has  already  been  shown  with 
Figures  5.22  to  5.25,  and  for  both  sets  of  tubes,  the 
optimum  fin  spacing  was  determined  to  be  1.5  mm.   Figure 
5.30  shows  the  enhancement  ratio  for  the  set  of  tubes  with 
2.0-mm  fin  height.   It  is  shown  that  the  maximum  enhancement 
is  about  5  and  6.5  for  the  vacuum  and  atmospheric  runs, 
respectively.   Figure  5.31  is  a  cross  plot  showing  the 
enhancement  ratio  versus  the  fin  spacing  for  both  pressure 
conditions,  with  the  fin  height  as  a  parameter.   It  is  shown 
that  the  set  of  tubes  with  2  mm  fin  height  has  better 
performance  than  the  set  of  tubes  with  1.0  mm  fin  height. 
The  reasoning  is  two  fold: 
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1.  Surface  area  increases  with  increasing  fin  height 
(see  Figure  5.30)  . 

2.  The  retention  angle  decreases  and  hence  heat-transfer 
performance  increases,  with  increasing  fin  height  as  pre- 
dicted by  equation  (1.12) .   This  was  also  visually  veri- 
fied during  this  thesis  effort. 

Examining  carefully  the  above  mentioned  Figure  5.31, 
it  can  be  seen  that  the  performance  sequence  is  not  the  same 
for  both  sets  of  tubes.   The  tube  with  1.0-mm  fin  spacing 
is  the  second  in  performance  among  the  set  of  tubes  with 
1.0-mm  height.   However,  for  the  set  of  tubes  with  2.0  mm 
height,  the  tube  with  1.0-mm  fin  spacing  is  the  third  in 
performance  while  the  tube  of  2.0-mm  fin  spacing  becomes 
second.   This  observation  holds  for  both  pressure  conditions 
This  means  that  as  the  fin  height  increases,  there  is  a 
tendency  for  the  optimum  spacing  to  move  to  a  higher  value. 
This  tendency  is  mainly  due  to  the  increased  amount  of 
condensate  which  is  retained  between  the  fins  at  the  top 
part  of  the  tube,  caused  by  the  increased  condensate 
generation.   This  fact  was  verified  visually  during  the 
experiments,  and  Figure  5.3  2  shows  a  schematic  of  the  con- 
densate appearance  on  those  two  tubes  which  differ  only 
by  the  fin  height.   It  is  shown  that  the  interfin  space  at 
the  top  part  of  the  tube  with  1.0-mm  fin  height  'holds'  less 
condensate  than  the  tube  with  2.0-mm  fin  height.   This  obser- 
vation may  reveal  that  the  mechanism  of  thinning  the  film 
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diminishes  when  the  fin  height  increases.   The  pressure 
gradient  in  the  condensate  film  between  points  A  and  B 
(see  Figure  1.1)  can  be  approximated  as: 


P-)  =   ^  (5.13) 

dy  I  ab         e 


Thus,  the  pressure  gradient  between  points  A  and  B  decreases 
with  increasing  fin  height  (e) ;  therefore,  it  should  be 
expected  that  the  tubes  with  2.0-mm  fin  height  would  provide 
a  less  effective  thinning  effect.   To  prove  this  statement, 
one  more  cross  plot  was  prepared  by  plotting  the  normalized 
enhancement  ratio  versus  the  fin  spacing  (i.e.,  separate 
the  area  addition  effect).   Figure  5.33  shows  clearly  that 
the  set  of  tubes  with  2.0-mm  fin  height  is  not  improving  the 
condensation  enhancement  mechanisms  as  effectively  as  "the 
set  of  tubes  with  1.0-mm  fin  height. 

It  is  worth  mentioning  that  as  the  fin  height 
increases,  the  portion  of  the  fin  that  extends  above  the  top 
flooded  part  of  the  tube  increases  (see  Figure  5.32  and 
compare  the  extending  fin  part  for  the  tubes  with  1.0  and 
2.0  mm  fin  height);  therefore,  this  part  of  the  fin  presents 
the  highest  heat-transfer  performance  since  the  condensate 
layer  around  the  fin  is  extremely  small.   Of  course,  the 
performance  of  the  tip  is  limited  by  the  fin  efficiency 
which  decreases  with  increasing  fin  height;  therefore,  the 
set  of  tubes  with  2.0-mm  fin  height  has  better  performance 
mainly  because  of  the  excessive  area  addition. 
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Up  to  now,  two  ways  to  increase  the  external  surface 
area  of  a  tube  were  mentioned;  i.e.,  either  by  reducing  the 
fin  thickness  or  by  increasing  the  fin  height.   In  the 
previous  section  it  was  shown  that  by  changing  the  fin 
thickness  there  was  not  much  enhancement,  the  reason  being 
that  the  fin  thickness  does  not  improve  the  flooding  problem 
and  that  the  area  addition  happens  in  the  regimes  of  the 
tube  that  are  flooded.   On  the  other  hand,  by  increasing  the 
fin  height  there  is  the  benefit  of  area  addition  away  from 
the  flooded  part  of  the  tube  resulting  in  a  better  heat- 
transfer  performance . 

5 .   Significance  of  Tube  Diameter 

Yau  et  al.  [Ref.  3]   took  data  for  filmwise  conden- 
sation on  copper  finned  tubes  under  near-atmospheric  pres- 
sure.  The  fin  shape  was  rectangular  and  the  fins  were  0.5 
mm  in  thickness,  1.0  mm  in  height,  and  the  fin  spacings  were 
from  0.5  up  to  20.5  mm  (13  different  fin  spacings).   The 
diameter  of  the  tubes  (at  the  root  of  fins)  was  12.7  mm  and 
the  vapor  velocity  was  1.1  m/s . 

Figure  5.34  shows  the  enhancement  ratio  versus  fin 
spacing  for  the  data  taken  by  Yau  et  al.  and  the  data  taken 
during  this  thesis  effort.   The  only  difference  between  the 
two  sets  of  tubes  was  the  root  diameter,  which  was  12.7  and 
19  mm  for  Yau  et  al .  and  the  present  work,  respectively. 
It  can  be  seen  that  there  is  a  similarity  in  the  enhancement 
ratio  curves  and  for  both  sets  of  tubes  the  best  fin  spacing 
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found  was  1.5  mm.   However,  there  is  a  significant  differ- 
ence in  the  enhancement  ratio,  which  can  be  explained  by 
the  following  reasons: 

1.  The  retention  angle  increases  as  the  tube  diameter 
decreases,  as  it  is  predicted  by  equation  (1.12) ,  resulting 
in  a  lower  condensing  coefficient. 

2.  The  smooth-tube  performance  is  predicted  by  the 
Nusselt  equation,  and  it  gives  that  a  tube  with  smaller 
diameter  will  have  higher  condensing  coefficient  (recall 
that  the  tube  diameter  is  in  the  denominator) . 

Therefore,  when  the  enhancement  ratio  is  computed, 
the  tube  with  smaller  diameter  will  have  the  numerator  small 
and  the  denominator  large,  resulting  in  a  smaller  ratio 
than  that  which  was  obtained  during  this  thesis  effort. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

1.  For  the  best  fin  spacing  of  1.5  mm,  a  fin  thickness 
of  1.0  mm  and  a  fin  height  of  2.0  mm,  heat-transfer  enhance- 
ments around  4  to  5  were  realized  for  the  vacuum  runs  and 
around  6  for  the  atmospheric  runs,  while  the  area  increment 
was  only  about  3  times  the  smooth-tube  area.   The  additional 
enhancement  is  due  to  the  thinning  effect;  therefore,  any 
theoretical  model  that  neglects  this  effect  may  have  severe 
deviations  from  the  experimental  data. 

2.  The  1.5  mm  fin  spacing  appears  to  be  the  optimum 
among  the  fin  spacings  tested  (0.5,  1.0,  1.5,  2.0  and  4.0  mm) 
for  both  vacuum  and  near-atmospheric  pressure  conditions. 

3.  The  condensation  enhancement  does  not  strongly  depend 
on  the  fin  thickness.   However,  the  tubes  with  fin  thickness 
around  0.75  to  1.0  mm  performed  better  than  the  tubes  with 
0.5  and  1.5  mm  fin  thicknesses. 

4.  As  the  fin  height  increased  from  1.0  to  2.0  mm, 

the  condensing  coefficient  increased  only  by  about  20  percent, 
despite  an  area  increase  of  about  50  percent.   This  lower- 
than-expected  performance  is  believed  to  be  owing  to  the 
poorer  thinning  effect  of  the  condensate  film  with  increas- 
ing fin  height. 

5.  As  the  tube  diameter  increases,  the  retention  angle 
decreases  resulting  in  a  better  heat-transfer  performance; 
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therefore,  the  tubes  with  a  root  diameter  of  19.0  mm  showed 
a  superior  enhancement  over  the  tubes  with  a  12.7  mm  root 
diameter. 

6.  The  consistency  of  equation  (1.12)  that  predicts  the 
retention  angle  was  visually  verified  throughout  all  the 
experiments.   Certainly,  it  must  be  a  very  useful  equation 
to  predict  the  condensate  retention  angle  when  modelling 
steam  condensation  on  finned  tubes. 

7.  The  leading  coefficient  of  the  Sieder-Tate-type 
equation  used  to  determine  the  water-side  resistance  with 
the  instrumented  tube  was  0.0635,  while  the  modified  Wilson 
plot  technique  gave  the  same  coefficient  as  0.071. 

8.  The  test  apparatus  successfully  operates  under  any 
pressure  condition.   An  extremely  tight  apparatus  is 
achieved  resulting  in  minimization  of  the  non-condensing 
gas  effect. 

9.  A  tube-cleaning  and  preparation  procedure  has  been 
established  to  ensure  filmwise  condensation  on  copper  tubes. 

B.   RECOMMENDATIONS 

1.  Calibrate  the  existing  pressure  transducer  in  order 
to  replace  the  U-tube,  mercury-in-glass  manometer.   This 
will  make  each  run  faster  and  more  automatic  (less  data  to 
be  entered  manually) . 

2.  Install  the  relief  valve  on  the  test  condenser. 

3.  Take  data  with  different  shaped  fins  (at  the  tip  as 
well  as  at  the  root  of  the  fin)  to  investigate  the  effect 
of  fin  geometry. 

120 


4.  Take  data  with  tubes  that  are  fabricated  with 
different  thermal  conductivity  materials  to  investigate  the 
effect  of  the  material  conductivity  on  the  condensing 
heat-transfer  coefficient. 

5.  Attach  drainage  strips  on  the  existing  tubes  to 
investigate  the  effect  of  the  drainage  strips  on  steam 
condensation  enhancement. 

6.  Take  data  to  study  the  effect  of  vapor  shear  on 
finned  tube  enhancement.   For  the  vacuum  runs,  it  is  recom- 
mended that  vapor  velocities  of  1.0,  2.0,  3.0  and  4.0  m/s 
be  used  and  for  the  atmospheric  runs  the  vapor  velocities 
of  0.7  and  1.0  m/s.   Higher  vapor  velocities  would  be 
desirable  if  the  apparatus  could  be  modified. 

7.  Modify  the  existing  apparatus  to  study  the  conden- 
sate inundation  effect  and  take  data  to  investigate  if  the 
enhancement  and  the  optimum  fin  spacing  are  affected  by 
inundation. 

8.  Take  data  on  finned  tubes  that  are  coated  with  a 
dropwise  promoter.   Since  the  presence  of  the  fins  will 
result  in  a  more  effective  bonding  of  the  promoter,  this 
combination  may  result  in  superior  enhancement  and  promoter 
endurance.   Also,  the  condensate  may  not  be  retained  on 
dropwise-coated  finned  tubes  with  sufficient  fin  spacing. 

9.  Modify  equation  (1.14)  of  Owen  et  al.  [Ref.  8] 

to  include  the  surface-tension  effect  on  the  unflooded  part 
of  the  tube.  By  including  this  effect,  the  Owen  prediction 
for  the  steam  data  will  become  less  conservative. 
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APPENDIX  A 
PREPARATION  OF  CONDENSER  TUBES 

Before  every  run,  the  tube  was  prepared  by  following 
the  steps  below.   The  operator  must  wear  gloves  to  avoid 
contamination  of  the  tube  as  well  as  to  avoid  skin  irritation 
that  may  be  caused  by  the  cleaning  solution. 

1.  Prepare  a  solution  of  equal  parts  of  sodium  hydroxide 
and  methyl  alcohol.   Heat  the  solution  until  it  becomes 
watery.   Using  a  toothbrush,  apply  a  coating  of  this  mixture 
on  the  tube. 

2.  Place  the  coated  tube  in  a  steam  bath  for  about  one 
hour. 

3.  Rinse  the  tube  with  tap  water  to  avoid  contamination 
of  the  apparatus  with  the  solution. 

4.  Clean  the  external  surface  of  the  tube  with  a  brush 
using  a  weak  solution  of  tap  water  and  a  detergent  (Spark- 
leen) .   This  will  ensure  that  all  the  sodium  hydroxide 
deposits  are  removed  from  the  corners  of  the  fins. 

5.  Using  a  brush,  clean  the  inside  surface  of  the  tube. 
It  is  essential  to  remove  any  fouling. 

6.  Install  the  tube  as  soon  as  possible  to  avoid  any 
system  contamination  from  the  environment. 
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APPENDIX  B 
SYSTEM  START-UP  AND  SHUT-DOWN  PROCEDURES 

To  start  the  system,  the  following  procedure  is 
recommended : 

1.  Ensure  that  the  boiler  water  level  in  the  boiler  is 
at  least  six  inches  above  the  upper  level  of  the  heaters. 

2.  Close  the  venting  valve  and  open  the  valve  on  the 
purge  line  (see  Figure  2.1). 

3.  Ensure  that  the  vacuum  pump  is  filled  with  oil  at  the 
right  level. 

4.  Turn  on  the  vacuum  pump. 

5.  If  the  run  is  for  low  pressure  then  turn  on  the  boiler 
power  when  the  system  pressure  approximates  the  desired 
pressure.   Avoid  boiling  if  the  system  pressure  is  greater 
by  50  mmHg  than  the  desired  pressure.   Boiling  at  an  ele- 
vated pressure  causes  the  apparatus  walls  to  reach  a 
temperature  higher  than  the  desired  saturation  temperature. 
If  this  happens,  a  considerable  time  interval  (as  much  as  30 
minutes  for  a  temperature  difference  of  10  K)  will  be  needed 
to  avoid  superheated  steam.   Note  that  when  the  system  walls 
are  at  an  elevated  temperature,  the  steam  gains  heat,  thus 
causing  superheating.   The  power  input  is  controlled  by  a 
potentiometer  until  the  voltmeter  reads  90  units  (note  that 
each  unit  represents  2.2  Volts).   This  setting  yields  a  2.0 
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m/s  vapor  velocity  at  a  system  pressure  of  85  mmHg .   If  the 
run  is  an  atmospheric  one,  then,  turn  on  the  power  right 
after  the  vacuum  pump  is  energized.   For  a  vapor  velocity 
of  1.0  m/s,  the  voltmeter  should  be  set  at  its  maximum 
which  is  170  units  or  (170.0  x2.2)  Volts. 

6.  Open  the  valve  to  fill  the  tap  water  sump. 

7.  Circulate  tap  water  through  the  auxiliary  condenser. 

8.  Check  that  the  water  sump  is  filled  and  then  turn  on 
the  pumps  that  circulate  water  through  the  test  bue. 

9.  When  the  steam  reaches  the  tube  (as  can  be  seen  from 
the  viewport) ,  adjust  the  water  flowrate  through  the  auxili- 
ary condenser  in  order  to  achieve  the  desired  pressure. 

10.   Energize  the  data  acquisition  system,  the  computer, 
the  printer  and  insert  in  the  driver  of  the  computer  the 
disk  which  contains  the  code  "DRP4."   Finally,  load  "DRP4" 
and  follow  the  interactive  program  to  gather  and  analyze 
data . 

To  secure  the  system: 

1.  Turn  off  the  input  power  to  the  boiler. 

2.  Secure  the  water  supply  to  the  auxiliary  condenser. 

3.  Open  the  venting  valve  to  relieve  the  pressure  in  the 
apparatus . 

4.  Turn  off  the  printer,  the  data  acquisition  system 
and  the  computer. 

5.  Secure  the  pumps  that  circulate  tap  water  through 
the  test  tube. 
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APPENDIX  C 
RAW  DATA  OF  UNINSTRUMENTED  AND  INSTRUMENTED  TUBE 

The  following  two  tables  contain  the  raw  data  taken  for 
the  uninstrumented  and  instrumented  tubes  under  vacuum  and 
near-atmospheric  conditions. 

TABLE  V 
Raw  Data  of  Uninstrumented  Smooth  Tube 


File 

Name 

:   WSIV14 

Pressure 

Condition 

:   Vacuum  (85  mmHg) 

Vapor 

Velocity 

:   1 

0  m/s 

Vw 

Tin 

Tout 

Ts 

m/s 

°C 

°C 

°C 

4.47 

18.58 

19.37 

49.45 

3.93 

18.63 

19.48 

49.59 

3.38 

18.70 

19.65 

49.46 

3.10 

18.71 

19.72 

49.76 

2.71 

18.81 

19.92 

49.41 

1.99 

18.99 

20.34 

49.84 

1.48 

19.19 

20.87 

49.49 

1.15 

19.46 

21.46 

49.53 

File 

Name 

-   WSIA16 

Pressure 

Condition 

:   Atmospheric 

Vapor 

Velocity 

:   0. 

8  m/s 

Vw 

Tin 

Tout 

Ts 

m/s 

°C 

°C 

°C 

4.49 

18.84 

20.65 

99.81 

3.93 

18.91 

20.94 

99.96 

3.38 

18.99 

21.31 

100.11 

3.10 

19.00 

21.48 

99.93 

2.71 

19.11 

21.85 

100.31 

1.99 

19.32 

22.84 

100.26 

1.48 

19.56 

23.92 

99.98 

1.15 

19.91 

25.20 

99.97 
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APPENDIX  D 
RAW  DATA  OF  TUBES  TESTED 

The  following  pages  contain  tables  with  the  raw  data 
taken  for  all  tubes  tested  under  vacuum  and  near-atmospheric 
conditions . 
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APPENDIX  E 
ERROR  ANALYSIS 

The  numerical  data  taken  during  an  experiment  are  used 
together  with  theoretical  formulations  in  order  to  obtain 
results  about  the  objective  of  the  experiment.   Any  measure- 
ment has  an  uncertainty  which  is  a  function  of  the  device 
accuracy  as  well  as  the  experience  of  the  operator. 
Therefore,  the  final  result  may  be  severely  distorted  due 
to  the  error  propagation  during  the  computational  steps. 
In  cases  where  the  final  result  has  a  large  uncertainty, 
it  may  be  unwise  to  draw  any  conclusions. 

The  probable  uncertainty  in  a  computation  can  be 
determined  from  the  following  equation,  as  taken  from 
[Ref .  19]  : 


1/2 
r  ,  3  R     .  2  ,  ,  8  R     »  2        ,  ,  9  R     \  2  -, 

R         dX.   1        3x~   2  3x„   n 


n 


(E.l) 


where 


R   is  the  result 

x-.  ,x»  ,  .  .  .  ,x   are  the  measured  (independent) 
1   2  .   '  n  c 

variables 

w, ,w_,...,w   are  the  uncertainties  in  the 
measured  (independent)  variables. 


The  above  equation,  together  with  the  uncertainty  of 
each  measurement,  was  used  to  determine  the  possible  error 
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of  the  outside  heat-transfer  coefficient.   As  discussed  in 
Chapter  V,  the  outside  resistance  was  found  by  subtracting 
the  water-side  and  tube-wall  resistances  from  the  overall 
resistance.   This  means  that  the  error  in  the  condensing 
coefficient  is  mainly  a  function  of  the  overall  and  the 
water-side  resistances.   In  such  cases,  if  the  water-side 
resistance  has  a  significant  error  it  will  be  unavoidably 
reflected  in  the  condensing  heat-transfer  coefficient. 
Further,  if  the  value  of  the  water-side  resistance  is  much 
larger  than  the  condensing  resistance,  then  the  error  that 
will  appear  in  the  condensing  coefficient  will  be  highly 
amplified.   The  steam-side  resistance  is  generally  smaller, 
especially  on  finned  tubes,  then  the  water-side  resistance 
(recall  that  the  condensation  process  is  a  phase  change 
phenomenon  where  the  heat-transfer  process  is  very  effec- 
tive) .   Therefore,  to  'protect'  the  final  result  of  the 
condensing  coefficient  from  an  excessive  error,  the  water- 
side and  condensing  resistances  should  be  as  nearly  balanced 
as  possible.   The  water-side  resistance  can  be  reduced 
either  by  increasing  the  coolant  velocity  or  by  internal 
enhancement.   Both  ways  promote  turbulence  and  reduce  the 
boundary  layer  thickness  on  the  inside  wall  surface  (recall 
that  the  coiled  insert  was  used  exactly  for  this  purpose) . 
However,  since  the  water-side  temperature  rise  decreases 
with  increase  in  water  velocity,  the  water  velocity  must  not 
be  too  excessive.   Since  the  water-side  temperature  rise  is 
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the  most  crucial  measurement,  a  minimum  value  must  be 
assigned  to  minimize  experimental  errors.   For  this  reason, 
the  internal  enhancement  is  preferable  than  the  increase 
in  water  velocity  to  obtain  a  lower  thermal  resistance  on 
the  waterside. 

Based  on  the  above  discussion,  it  should  be  expected 
that  as  the  cooling  water  flow  rate  decreases  (i.e.,  less 
turbulence) ,  the  error  in  the  outside  heat-transfer  coeffi- 
cient should  increase.   Figures  E.l  and  E.2  show  the  error 
bands  for  three  cooling  water  flow  rates  for  the  tube  that 
showed  the  best  heat- transfer  performance  (1.0  mm  fin  thick- 
ness, 1.5  mm  fin  spacing  and  2.0  mm  fin  height)  for  both 
pressure  conditions.   It  is  shown  that  indeed  the  error 
gets  larger  as  the  cooling  water  flow  rate  (i.e.,  the  heat 
flux)  decreases.   Also  notice  that  the  error  bands  shown 
with  the  atmospheric  run  are  smaller  than  the  corresponding 
value  for  the  vacuum  runs.   This  improvement  can  be  ex- 
plained by  the  increased  water-side  temperature  rise. 
Note  that  the  relative  error  decreases  as  the  water-side 
temperature  rise  increases. 

Based  on  the  foregoing  discussion,  it  is  clear  why  the 
reference  heat  fluxes  for  computing  the  enhancement  ratio 

of  the  cross  plots  in  Chapter  V  were  selected  as  high  as 

2 

possible  (0.35  and  1.0  MW/m   for  vacuum  and  atmospheric 

runs,  respectively). 

The  uncertainty  of  the  major  variables  that  govern  the 
condensation  process  (cooling  water  velocity,  Reynolds 
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number,  heat  flux,  LMTD,  water-side  and  outside  heat- 
transfer  coefficients)  are  given  as  percentages  in  the 
following  computer  outputs  of  the  program  "ERROR. "   At 
the  end  of  this  appendix,  the  listing  of  the  program 
"ERROR"  is  provided  for  the  Error  Analysis. 
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DATA  FOR  THE  ERROR  ANALYSIS: 

File  Name:  F25E2V146 

Pressure  Condition:   VACUUM  11  kPa 

Cooling  Water  Temperature  In   =  22.72  (Deg  C) 

Cooling  Water  Temperature  Out  =  24.07  (Deg  C) 

Steam  Temperature             -  48.35  (Deg  O 

Water  Flou  Rate  <%)           =  80.00 

Water  Velocity                =  4.39  (m/s> 


ERROR  ANALYS 


TC 


% 

Error 

in 

Mass  Floy 

Rate 

= 

1  .03 

% 

Error 

m 

Reynolds  1 

dumber 

= 

1  .24 

Error 

in 

Heat  Flux 

= 

0.79 

X 

Error 

in 

LMTD 

= 

2.22 

V. 

Error 

m 

Uo 

= 

1 

Error- 

in 

Hi 

= 

1  .32 

V 
/> 

Error 

in 

Ho 

= 

3 .  36 

4- 


DATA  FOR  THE  ERROR  ANALYSIS: 

File  Name:  F25E2V146 

Pressure  Condition:   VACUUM  M  kPa 
Cooling  Water  Temperature  In 
Cooling  Water  Temperature  Ou 
Steam  Temperature 
Water  Flou  Rate  <  7, ) 
Water  Ve loc i  ty 


ERROR  ANALYSIS: 

%  Error  in  Mass  Flou  Rate  =  1,80 

7-.    Error  in  Reynolds  Number  =  1  .93 

'A    Error  in  Heat  Flux  -  0.73 

V,   Error  m  LMTD  -  1  .S3 

%  Error  in  U o  ! . 8  I 

V,   Error  in  Hi  -  1  .78 

%  Error  in  Ho  =  '0,53 


9 ''-"! 

.02 

(Deg 

C) 

24 

.79 

(Deg 

C) 

48 

,20 

(  Deg 

C) 

45 

.00 

2 

,51 

<  m  /  s 

> 
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DATA  FOR  THE  ERROR  ANALYSIS 


r  i 


e  Name:  F25E2VI46 

Pressure  Condition:   VACUUM  11 
Cooling  Water  Temperature  In 
Cooling  Water  Temperature  Out  : 
Steam  Temperature 
Water  ^lou  Rate  <%) 
Water  Velocity 


Pa 

23 

.53 

( Deg 

C  > 

26 

.  !  1 

(Deg 

C) 

48 

53 

(Deg 

C) 

20 

.00 

1 

16 

(m/s) 

ERROR  ANALYSIS 


Error 

in 

Mass  Flow  Rate 

7 

Error 

in 

Reynolds  Number 

% 

Error 

m 

Heat  Flux 

7. 

Error 

in 

LMTD 

v 

Error 

in 

Uo 

7 

Error 

in 

Hi 

7 

Error 

i.  n 

Ho 

.39 

3 

.96 

« 

.30 

1 

.  10 

1 

.70 

29 

=  36.25 


DATA  FOR  THE  ERROR  ANALYSIS: 


F25E2A128 


rile  Name : 
Pressure  Condition:   ATMOSPHERIC 
Cooling  Water  Temperature  In 
Cooling  Water  Temperature  Out  = 


bream 


'mpe rat ure 


Water  Flow  Rate 
Water  Veloci  ty 


( X ) 


to 

'<Pa 

22 

90  (Deg  C> 

26 

96  (Deg  C> 

99 

37  (Deg  O 

80 

00 

A 

39  (m/s) 

CO 


RROR  ANALYSIS: 


X  Error 
X  Error 
X  Error 
7.  Error 
V,  Error 
7  E  r  r  o  r 
%  E  r  r  o  r 


in  Mass  Floy  Rate 
in  Reynolds  Number 
Heat  Flux 


m 
in 
in 
in 
in 


LMTD 
Uo 
Hi 
Ho 


-42 

71 
S3 
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DATA  FOR  THE  ERROR  ANALYSIS: 

File  Name:  F25E2AI2S 

Pressure  Condition:   ATMOSPHERIC  101  kPa 

Cooling  Water  Temper at are  In 

Cooling  Water  Temperature  Out  = 

Steam  Temperature 

Water  Flou  Rate  (%) 

Water  Velocity 


ERROR  ANALYSIS: 

7.   Error  in  Mass  Floy  Rate  =  I  .30 

V.    Error  in  Reynolds  Number  =  1  .95 

%  Error  in  Heat  Flux  =  0.61 

V.  Error  in  LMTD  -  0.51 


OO 

.27 

(Deg   C) 

28 

.37 

(Deg   C) 

99 

.59 

(Deg    C) 

45 

.00 

*-« 

,51 

(  m  /  s  ) 

30 


%   Error  in  Uo  -      U 

%  Error  in  Hi  1.73 

V,    Error  in  Ho  =   6.34 


DATA  FOR  THE  ERROR  ANALYSIS: 

File  Name:  F25E2AI28 

Pressure  Condition:   ATMOSPHERIC  !01  kPa 

Cooling  Water  Temperature  In   =  23.96  (Deg  O 

Cooling  Water  Temperature  Out  =  32.30  (Deg  C) 

Steam  Temperature             =  99.43  (Deg  C> 

Water  Flou  Rate  <%>           =  20.0  0 

Water  Velocity                =  1.16  (m/s> 


ERROR  ANALYSIS: 


Er 


ror  in  Mass  Flou  Rate 

y  Error  in  Reynolds  Number 

7,    Error  in  Heat  Flux 

y  Error  in  LMTD 

7.   Error  in  Uo 

K  Error  in  Hi 

X  Error  in  Ho  =    27  ,\)\ 


o 

.30 

"1 

,97 

1 

,25 

0 

,34 

1 
I 

.30 

"3 
O 

.30 
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! 000 !  FILE  NAME  :  ERROR 

1010!  REVISED  :  JUNE  27.  1984 

1020! 

1030  PRINT  USING  "1  OX .""DATA  FOR  THE  ERROR  ANALYSIS: 

1040  PRINT 

1050  BEEP 

10G0  INPUT  "ENTER  FILE  NAME".FileS 

1070  PRINT  USING  "1  OX  .  ""Fi  le  Name  :  "" . 1 2A" ;Fi le£ 

1080  BEEP 

!090  INPUT  "ENTER  PRESSURE  CONDITION  <  1  -VACUUM. 2-ATMOSPHERIC)  ".Pre 

1100  IF  Prc=1  THEN 

1110  BEEP 

1120  PRINT  USING  "1  OX  .  ""Pressure  Condition:   VACUUM  11  IcPa 

1130  ELSE 

I  1 40*  PRINT  USING  "TOX  .""Pressure  Condition:   ATMOSPHERIC  101  kPa 

1150  END  IF 

II  GO  BEEP 

1170  INPUT  "ENTER  COOLING  WATER  TEMPERATURE  IN". Tel 

1180  PRINT  USING  "1  OX ,""Cool ing  Water  Temperature  In   =  "" .3D.2D ,""  (Deg  CV" 

Tci 

1190  BEEP 

1200  INPUT  "ENTER  COOLING  WATER  TEMPERATURE  DIFFERENCE" ,Dtu» 

1210  Tco=Tci+Dtu 

1220  PRINT  USING  "TOX .""Cool ing  Water  Temperature  Out  =  "".3D. 2D,""  (Deg  C)"' 

Tco 

1230  BEEP 

1240  INPUT  "ENTER  STEAM  TEMPERATURE" . Ts 

1250  PRINT  USING  "10X. ""Steam  Temperature  =  "",3D.2D,""  (Deg  C)"' 

Ts 

1260  BEEP 

1270  INPUT  "ENTER  COOLING  WATER  FLOW  RATE  7."  ,Fn 

1280  PRINT  USING  "1  OX  ,  ""Water  Flow  Rate  (*/.)  =  ""  ,  3D  .  2D"  ;Fm 

1290  T=<Tci+Tco)/2   !  FILM  TEMPERATURE 

1300!  UNCERTAINTY  IN  THE  COOLING  WATER 

1310  Dfm-.IO  !  ERROR  IN  FLOW  METER  READING 

1320  Drho=1   !  ERROR  IN  WATER  DENSITY 

1330  Dmf=5.00049E-3+6.9861937E-3*Dfm  !  ERROR  IN  MASS  FLOW  RATE 

1340  Rho=FNRho(T>  !  WATER  DENSITY 

1350  Mf=1 ,04805E-2+6.80932E-3*Fm   !  MASS  FLOW  RATE  OF  COOLING  WATER 

13G0!  CORRECT  MF  FOR  THE  TEMPERATURE  EFFECT 

1370  Mf=Mf»( 1 .0365-1 .96644E-3*Tci+5 .252E-6*Tci *2>/ .995434 

1380  Di= .0127   !  TUBE  INSIDE  DIAMETER 

1390  Ai-<PI»Di'2>/4  !  TUBE  INSIDE  CROSS  SECTION  AREA 

1400  Ddi=. 000025 

1410  Dai=<PI*Ddi  2)/4  !  ERROR  OF  INSIDE  TUBE  CROSS  AREA 

1420!  COMPUTE  THE  WATER  VELOCITY 

1430  Vu=Mf/(Rho*Ai)  !  WATER  VELOCITY 

1440"  PRINT  USING  "10X, ""Water  Velocity  =    "".Z.DD.""  (m/s) 

W 

1 450!  CORRECT  OUTLET  WATER  TEMP.  FOR  THE  MIXING  CHAMBER  EFFECT 

1460  Tco=Tco-.004*VuT2  !  COOLING  WATER  OUTLET  TEMP. 

1470  T=(Tci+Tco)*.5  !  FILM  TEMPERATURE 

1480!  COMPUTE  THE  ERROR  IN  WATER  VELOCITY 

1490  Dvw-Vw*SQR<<Dmf/Mf ) 2+< Drho/Rho) '2+<Dai/Ai) "2> 

1500!  UNCERTAINTY  IN  THE  REYNOLDS  NUMBER 

1510  Mu-FNMui(T)  !  WATER  VISCOSITY 

1520  Dmw=6.E-6  !  ERROR  OF  WATER  VISCOSITY 

1530  Re=(Rho*Vu*Di)/Mu 

1540  Dre=Re*SGR( (Drho/Rho)  "2+  <  Dvui/Vw)    2  +  (Ddi/Di  )  '2  +  (Dmu/Mu)    2) 
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1550!  UNCERTAINTY  IN  THE  HEAT  TRANSFERED 

1560  Cpu=FNCpu<T) 

1570  Q=Mf*< Tco-Tci)*Cpu 

1580  Dtco=.02 

1590  Dtci=.02 

1600  Dcpu=8 

1610  Dq=Q*SQR(<Dmf/Mf >  2+< (Dtco/ < Tco-Tc i > > >  -2  +  < (Dtci/( Tco-Tci ) ) ) "2  +  < Dcpu/Cpui)  '2 

) 

1620!  UNCERTAINTY  IN  THE  HEAT  FLUX 

1630  Dl= .0005  !  ERROR  IN  TUBE  LENGHT 

1640  Do-. 01905  !  SMOOTH  TUBE  OUTER  DIAMETER 

1650  Ddo=. 000025 

1660  L=. 13335  !  FINNED  TUBE  LENGTH 

1670  Qp=0/(PI»Do*L)  !  HEAT  FLUX 

1680  Dqp«<Qp/PI>»SQR<<Dq/Q>*2  +  <Ddo/Dor2  +  <Dl/L>*2> 

1690  Dts=.02   !  ERROR  IN  TEMPERATURE  MEASUREMENT 

1 700  Lmtd=<Tco-Tci>/LOG(<Ts-Tci)/(Ts-Tco>> 

1710  Uo=Qp/Lmtd  !  OVERALL  HEAT  TRANSFER  COEF . 

1720  A1=Dts*(Tci-Tco)/((Ts-Tct)*(Ts-Tco)*L0G((Ts-Tci)/(Ts-Tco))) 

1730  A2=Dtci/((Ts-Tci)*L0G(<Ts-Tci)/(Ts-Tco))) 

17  40  A3=Dtco/<(Ts-Tco)*LGG(<Ts-Tci)/(Ts-Tco))) 

1750  Dlmtd=Lmtd*SQR(A1 "2+A2"2+A3 '2> 

1760  Du.o-Uo*SQR((Dqp/Qp)"2+<Dlmtd/Lmtd)  '2) 

1770  M=Mu 

1780  T1=<T+273.15)/273. 15 

1790  Ku=FNKw<T1) 

1800  Ci=.0635  !  SIEDER-TATE  COEFFICIENT 

1810  Ac=26.4   !  INTERSCEPT  FROM  SIEDER  PROGRAM 

1820  Kc=385  !  WALL  THERMAL  CONDUCTIVITY 

1830  L1-. 060325  !  LENGTH  OF  UNFINNED  LEFT  PART  OF  TUBE 

1840  L2=. 034925  !  LENGTH  OF  UNFINNED  RIGHT  PART  OF  TUBE 

1850  Dr=. 015875 

1360  Pr=Cpui*Miu/Kui 

1870  Mnu)=FNMuu<T> 

1380!  UNCERTAINTY   OF    INSIDE   HEAT-TRANSFER   COEFF. 

1890  Cf-1. 

1900  Hi-<Kw/Di>*<Ci»Re' .8»Pr  \333*Cf +Ac) 

1910  Dti=Q/(PI*Di*(L+L1*Fe1+L2*Fe2)*Hi) 

1320  Cfc=<MuWFNMuw<T+Dti)>  \  !4 

1930  IF   ABS((Cfc-Cf)/Cfc)>.01    THEN 

1940  Cf=(Cf+Cfc)*.5 

1950  GOTO    1900 

1960  END    IF 

1970  P1=PI»(Di+Do) 

1980  B1=<Do-Di)*PI*(Di+Do)*,5 

1990  M1-<Hi«P1/<Kc*B1 >>' .5 

2000  P2=PI*(Di+Dr) 

2010  B2=(Dr-Di)*PI*<Di+Dr)*.5 

2020  M2=(Hi*P2/(Kc*B2))  .5 

2030  Fe1=FNTanh(M1*L1 >/<M1»L1> 

2040  Fe2=FNTanh(M2»L2)/(M2*L2) 

2050  Dtc=0/(PI*Di*(L+L1*Fe1+L2*Fe2)*Hi) 

2060  IF  ABS((Dtc-Dti)/Dtc)>.01  THEN  1900 

2070  Dkw= . 001 0  !  ERROR  IN  WATER  THERMAL  CONDUCTIVITY 

2080  Dci=.0005  !  ERROR  IN  SIEDER-TATE  COEFFICIENT 

2090  DPr=.05  !  ERROR  IN  PRANDTL  NUMBER 

2100  Dcf=8.E-6 

2110  A4=.14*Dcf/Cf 

2120  Dhi-Hi*SQR<<Dlru/Ku>  -2  +  <Ddi/Di)    2+<  .3*Dre/Re>    2  +  <  .333*Dpr/Pr)    2+<Dci/Ci>  "2+ 

A4> 

2 '30!  UNCERTAINTY  OF  OUTSIDE  HEAT-TRANSFER  COEFF. 
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2140 

2150 

2160 

2170 

i))  )• 

2'80 

2190 

2200 

2210 

2220 

2230 

2240! 

2250 

2260! 

227,0 

2280! 

2290 

2300! 

2310 

2320! 

2330 

2340! 

2350 

2360! 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 


Dkc=1  !  ERROR  IN  TUBE  THERMAL  CONDUCTIVITY 
Ru=Do*L0G(Do/Di)/<2*Kc)  !  WALL  RESISTANCE 
Ho=1/( (1/Uo)-<Do*L/(Di*(L+L1*Fe1 +L2*Fe2)*Hi> )-Ru) 

Dru)=Ru*SGR<  (Ddo/Do)  '2  +  (Dlcc/Kc)  ~2  +  <  Ddo/(  Do*L0GCDo/Di ) )  >  '2+<Ddi/(Di*L0G<Do/D 
2) 
A5=!/Uo-Rui-(Do/(Di*Hi)) 
A6=Duo/<Uo"2*A5) 
A7=Dru/A5 

A8=<(Do/(Di*Hi))*<Dhi/Hi>>/A5 
PRINT 

Dho=Ho*SQR<A6~2+A7~2+A8~2> 
CALCULATE  THE  7.  ERROR  IN  Ho 
Prho=Dho*100/Ho 

CALCULATE  THE  7.  ERROR  IN  REYNOLDS  NUMBER 
Prre=Dre*100/Re 
CALCULATE  THE  7.  ERROR  IN 
Prmf  =  Dinf»100/Mf 
CALCULATE  THE  7.  ERROR  IN 

Prqp=Dqp*100/Qp 

CALCULATE  THE  7.  ERROR  IN 

Prlmtd=Dlmtd»100/Lmtd 

CALCULATE  THE  7.  ERROR  IN 

Pruo=DuLO*100/Uo 

CALCULATE  THE  7.  ERROR  IN 

Prhi=Dhi*100/Hi 

PRINT 

USING 


MASS  FLOW  RATE 

HEAT  TRANSFER 

LMTD 

OVERALL  HEAT  TRANSFER  COEF. 

INSIDE  HEAT  TRANSFER  COEFF. 


•10X,  ""ERROR  ANALYSIS: 


USING 
USING 
USING 
USING 
USING 
USING 
USING 


"10X, 
"10X 
"10X, 
"1  OX . 
"10X, 
"1  ox 

"I  OX. 


""7. 


Error 
Error 
Error 
Error 
Error 
Error 
Lrror 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

END 

DEF  FNMuu(T) 

A=247.8/<T+133.15> 

Muu=2.4E-5*10"~A 

RETURN  Muu 

FNEND 

DEF  FNTanh(X) 

P=EXP(X) 

0=EXP<-X> 

Tanh=<P-Q)/(P+Q) 

RETURN  Tanh 

FNEND 

DEF  FNKw<T1> 

Ku=-.92247+T1*<2.8395 

RETURN  Ku 

FNEND 

DEF  FNMw(T) 

A=247.8/(T+133.15) 

Mu=2.4E-5*10'A 

RETURN  Mu 

FNEND 

DEF  FNRho(T) 

Rho=999.52946+T*( .01269-W5 

RETURN  Rho 

FNEND 


in  Mass  Flow  Rate 

in  Reynolds  Number 

in  Heat  Flax 

in  LMTD 

in  Uo 

m  Hi 

in  Ho 


•\2D.2D.";Prmf 
".2D.2D.":Prre 
"".Z.2D.";Prap 
"*\Z.2D";Prlmtd 


i*  ft   •"* 


Z.2D.' 
.2D.2D. 
2D.  2D.' 


Pruo 
Prhi 
Prho 


T1*<1  .8007-TW .52577-.07344*T1>>) 


48251 3E-3-T*i .234147E-5)) 
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2730  DEF    FNCPw<T) 

2740  Cdu=<4.21 1 20358-T*<  2 . 26S2SE-3-T -<  4 . 4236 1 E-5+2 . 7 1 428E-7-T) ) >• I  000 

2750  RETURN    Cpu 

2750  FNEND 
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